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We consider methods to identify the classical ground state for an exchange-coupled Heisenberg antiferro- 
magnet on a non-Bravais lattice with interactions J 4j to several neighbor distances. Here we apply this to the 
unusual "octahedral" lattice in which spins sit on the edge midpoints of a simple cubic lattice. Our approach is 
informed by the eigenvectors of Jij with largest eigenvalues. We discovered two families of non-coplanar states: 
(i) two kinds of commensurate state with cubic symmetry, each having twelve sublattices with spins pointing 
in (1,1,0) directions in spin space (modulo a global rotation); (ii) varieties of incommensurate conic spiral. The 
latter family is addressed by projecting the three-dimensional lattice to a one-dimensional chain, with a basis of 
two (or more) sites per unit cell. 

PACS numbers: 75.25.-j, 75.30.Kz, 75.10.Hk, 75.40.Mg 



I. INTRODUCTION 

This paper concerns the classical ground state of the Hamil- 
tonian 



« = E - J * 



(i.i) 



where {s^} are unit vectors, and the couplings {Jij} have 
the symmetry of the lattice and may extend several neighbors 
away (being frustrated in the interesting cases). 

After an antiferromagnet's ordering pattern (or partial in- 
formation) is determined by neutron diffraction, the next ques- 
tion is which spin Hamiltonian(s) imply that order, if we admit 
interactions J 2 to second neighbors or J„ to further neighbors. 
The starting point for understanding ordered states is always 
the classical ground state(s). If the spins sit on a Bravais lat- 
tice (e.g. face-centered cubic), the solution is trivial due to a 
rigorous recipe, called the "Luttinger-Tisza" (LT) method (see 



Sec. II A below): the spins adopt (at most) a simple spiral - 
a coplanar state, meaning all spins point in the same plane 
of spin space JT] [2]. But if the spins form a lattice with a 
basis (more than one site per primitive cell), - e.g. kagome, 
diamond, pyrochlore, or half-garnet lattices - no mechanical 
recipe is known to discover the ground state. In these more 
complicated lattices, magnetic frustration (competing interac- 
tions) often induces complicated spin arrangements. 

Our aim has been to find a recipe for general lattices (albeit 
neither exhaustive nor rigorous) to discover the ground state 
spin pattern corresponding to a given set of exchange cou- 
plings Ji, to neighbors at successive distances. That is obvi- 
ously a prerequisite for solving the inverse problem (given the 
ordering patterns found by neutron diffraction, which com- 
bination^) of interactions can explain them?). Furthermore, 
after the whole phase diagram is mapped out, we can iden- 
tify the parameter sets leading to exceptionally degenerate or 
otherwise interesting states, so as to recognize which real or 
model systems might be close to realizing those special states. 

In this work, we focus on a narrower question: which pa- 
rameter combinations give a noncoplanar ground state, which 
could never happen in a Bravais lattice? We adopt the ex- 
change Hamiltonian \\.\\ , excluding single-site anisotropies 



and Dzyaloshinskii-Moriya couplings, which can trivially 
give non-coplanar ground states. 

We do not count cases where a non-coplanar ground state 
belongs to a degenerate family of states that also includes 
coplanar ground states. This happens trivially when two sub- 
lattices aren't coupled at all, or nontrivially when the inter- 
actions are constrained to cancel. In the latter cases thermal 
or quantum fluctuations usually break the degeneracy, favor- 
ing the collinear or coplanar states [3 , 4|. (A small amount of 
site-dilution or bond disorder can generate a uniform effective 
Hamiltonian that favors non-coplanar states (3]|5l, but here we 
only consider undisordered systems.) 



Motivations for noncoplanarity 

There are specific physical motivations to hunt for non- 
coplanar states. First, they point to possible realizations of 
chiral [6| spin liquids, such as are described within bosonic 
large-A r formalisms (as are hoped to approximate the behav- 
ior of frustrated magnets with s = 1/2). Such formalisms 
describe transitions from an ordered state to a quantum- 
disordered spin liquid; since there is no generic reason for 
a state to stop being chiral at the same time it loses spin or- 
der, a chiral ordered state presumably transitions into a chiral 
spin liquid. Hence, as a rule of thumb, a chiral spin liquid is 
feasible if and only if the classical ordered state (on the same 
lattice) is non-coplanar. 00. 

Secondly, spin non-coplanarity in metals (usually induced 
by an external magnetic field) allows the anomalous Hall ef- 
fect observed in pyrochlore and other magnets. Il9i fi"2"l This is 
ascribed to spin-orbit coupling and the Berry phases of hop- 
ping electrons (which are zero in the collinear or coplanar 
case). 

Thirdly, the symmetry-breaking of noncoplanar exchange- 
coupled magnetic states is labeled by an order parameter 
which is an 0(3) matrix, so the order-parameter manifold is 
disconnected. This permits a novel topological defect: the 
Z2 domain wall [13|, which is only possible in non-coplanar 
phases. 

Finally, there is current interest in "multiferroic" materials 
(i.e. those with cross couplings of electric and magnetic polar- 
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izations). For example, in the canonical multiferroics RMn0 3 
(where R=rare earth), frustrated exchange interactions induce 
a coplanar spiral, which in the presence of Dzyaloshinskii- 
Moriya anisotropic interactions carries an electric polariza- 
tion with it |[T4l4T8ll . If these spirals were asymmetric conic 
spirals, like our second class of ground states, there is gener- 
ically a net moment along the axis, which serves as a conve- 
nient "handle" to externally manipulate the orientation of the 
ground state (and thus control the multiferroic properties). 



A. The octahedral lattice 

Our spins sit on a rarely studied lattice we christen the "oc- 
tahedral lattice", consisting of the medial lattice (bond mid- 
points) of a simple cubic lattice, thus forming corner-sharing 
octahedra. Thus, each unit cell has a basis of three sites, 
forming what we call the x, y, and z sublattices (according 
to the direction of the bond they sit on). Each cubic ver- 
tex is surrounded by an octahedron of six sites, with nearest- 
neighbor bonds forming its edges; these octahedra share cor- 
ners, much as triangles or tetrahedra share corners in the well- 
known kagome and pyrochlore lattices. (Indeed, although the 
"checkerboard" lattice was introduced as a two-dimensional 
version of the pyrochlore lattice 03D, the octahedral lattice 
is the best three-dimensional generalization of the checker- 
board lattice.) This lattice was first studied as a frustrated 
Ising antiferromagnet ll20l |2"T1 More recently, it was used 
as a (simpler/pedagogical) toy model in papers aimed at the 
"Coulomb phase" of highly constrained spins on a pyrochlore 
lattice l22l l23l . It is one of the lattices constructed from the 
root lattices of Lie algebras. Il24l 

The octahedral sites are Wyckoff positions (and hence can- 
didates for a magnetic lattice) in most cubic space groups, 
so this is plausible to find in real materials, and a few are 
known. Most simply, it is realized by the transition metal 
sites in the CU3AU superstructure of the fee lattice ll20l (i.e. 
all but one of the four simple cubic sublattices). So far, the 
only example known in which the "Cu" lattice is magnetic 
seems to be Mn 3 Ge which is ferromagnetic [25] Another re- 
alization is in metallic perovskites such as Mn3SnN, though 
again the known materials are ferromagnetic ll26l . It would 
actually seem quite plausible to find realizations of our mod- 
els, which mostly have several exchange interactions Ji with 
competing signs, among metallic alloys: the RKKY interac- 
tion, expected between local moments in any metal, oscillates 
with distance inside a slowly decaying envelope. 

The octahedral lattice is closely related to the magnetic lat- 
tice found in the (mostly metallic) Ir 3 Ge 7 structures, including 
the strong-electron-interaction superconductor Mo3Sb7 ll27l 
[28). In that lattice, the simple-cubic lattice sites are sur- 
rounded by disjoint octahedra, i.e. a dimer of two magnetic 
ions decorates each bond of the simple cubic lattice. If this 
dimer were strongly coupled ferromagnetically, it would be a 
good approximation to treat it as a single spin, which is ex- 
actly the octahedral lattice. Instead, in MoaSby the dimers 
are antiferromagnetically coupled and, since Mo has spin 1/2, 
they form singlets l28l . If the spin length were longer, justify- 



ing classical treatment, we could convert to the ferromagnetic 
case simply by inverting the spin directions in every octahe- 
dron around an odd site of the cubic lattice, and changing the 
sign of all bonds coupling even sites with odd sites. Thus, 
much of the classical phase diagram for the MosSby lattice is 
related to that of the octahedral lattice. 

In this paper, we mainly consider four kinds of couplings, 
for separations out to the third neighbors: Ji for (1/2, 1/2, 0), 
J 2 or J' 2 for (1,0,0), J 3 for (1, 1/2, 1/2). Notice that cou- 
plings with the same displacement need not be equivalent by 
symmetry, since the site symmetry is just fourfold, less than 
cubic. (Our naming convention is to use the prime for the 
separation which requires more first-neighbor steps to tra- 
verse.) We also (less extensively) consider interactions J4 
or J4 for (1,1,0). To organize our exploration of this pa- 
rameter space, in analytic calculations we shall often assume 
J3, J4,, J4 <C J%, J2, J'2 (that suffices to give examples of most 
of the classes we found of noncoplanar ground states). 



B. Outline of paper and preview of results 



We begin (Sec. [TTJ by developing the techniques and con- 
cepts necessary to find the phase diagram as a function of 
the Ji/s and to discover non-coplanar ground states. We 



found ground states using three methods. The first (Sec. II A 
was Fourier analysis, known as the "Luttinger-Tisza" method, 
which can give a lower bound on the energy, but may not give 
a full picture of the ground state. The second (Sec. II B is an it- 



erative minimization algorithm, which numerically converges 
to a ground state; we introduce several diagnostic tools for un- 
derstanding the spin patterns produced by iterative minimiza- 



tion. The third method (Sec. II C 1 is the variational optimiza 



tion of idealized patterns displayed by iterative minimization. 

We then turn to our results, beginning with descriptions of 
the several classes of magnetic state we found for the octahe- 



dral lattice: various coplanar states (Sec. Ill I, the noncopla- 
nar, commensurate "cuboctabedral" spin states (Sec. |IV) , and 
a more generic group of noncoplanar, incommensurate "conic 
spirals" (Sec. [V]); in these the lattice breaks up into layers 
of spins with the same directions, each layer being rotated 
around the same (spin-space) axis relative to the layer be- 
low. Particularly noteworthy was a "double-twist" state we 
encountered, which is something like a conic spiral which 
also has a complex modulation in the transverse directions 



(Sec. VI 1. The plain stacked structures can be studied by 
mapping to one-dimensional ("chain") lattices, also with cou- 



plings to many neighbors, as worked through in Sec. V B 



From this we go on (Sec. VII 1 to quickly survey the phase 



diagrams we found, first for the cuboctahedral lattice, and 
then for the chain lattice (when treated as a lattice in its own 
right). In the conclusion, Sec. VIII we reflect on what our 
results might suggest for other frustrated lattices, such as the 
pyrochlore. 
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II. METHODS AND FRAMEWORK 

We employed several approaches to discover and under- 
stand ground states, for each given set of interactions (these 
are developed rest of this section - except for the use of map- 
pings, which we explain in Sec. [V] where it becomes natural 
to employ this technique). 

(a) . A Fourier analysis of the Hamiltonian (JTTTJ as 
a quadratic form with coefficients J^, the so-called 



"Luttinger-Tisza" method outlined in Sec. II A 



(b) 



. Iterative minimization, our main "exploratory" tech- 
nique. Starting from a random initial condition, we suc- 
cessively adjusted randomly chosen spins so as to re- 



duce the energy. (Sec. II B i. We then analyzed each re- 
sulting pattern with various diagnostics, as described in 
Sec. |IIB| and tagged the non-coplanar ones for further 
investigation. 

(c) . Variational optimization of the iterative minimization 
ground state. Finding a closed-form for the ground state 
introduces a number of free parameters (the most obvi- 
ous being a wave-vector). By allowing these parameters 
to vary from the values found with iterative minimiza- 
tion, we find a new, more rigorous, ground state. 

(d) . Mapping the (three-dimensional) problem to a simi- 
lar problem in a one-dimensional "chain" lattice with a 
basis of two sites. This is valid when the optimal (three 
dimensional) spin configuration is a stacking of layers, 
which we judged based on the results from approaches 
(a) and (b). The states on this simplified chain lattice 
may be found using approaches (a) and (b), or analyt- 
ically solved after parametrizing the state with a set of 
variational parameters. 



A. Spin states and Luttinger-Tisza modes 

The general theory of spin arrangements is reviewed in 
Refs. |2] and [29j The most fruitful approach to finding the 
ground states of the Hamiltonian dTTTT ) is to treat it as a 
quadratic form rewriting as 



H = - 



J2 J«^(k)§ a (-k) • ^(k), (2.1) 
where a and /? are sublattice indices; the explicit formulas 



k a,0=l 



for the cuboctahedral lattice case are given in Appendix A 1 
Eq. ( |A 1 [ >. Then we diagonalize this matrix, obtaining 



(2.2) 



k v=l 



Here m — 3 is the number of sites per primitive cell, and v is 
a band index; thus { J(kzy)} are the eigenvalues of Jy as an 
mN x mN matrix, with N being the number of cells (to be 
taken to infinity), and the wavevector k runs over the Brillouin 



zone. [In the Bravais lattice case, m — 1, the eigenvalue J(k) 
is simply the Fourier transform of J(r); for m > 1.] Also, 
s(k^) (complex-valued 3-vector) is the projection of the spin 
configuration onto the corresponding normalized eigenmode, 
A r_1 / 2 w l/s exp i(k ■ r). We shall call these the "Luttinger- 
Tisza" (LT) eigenvalues and modes [2, 30, 31 1; a mode with 
the most negative J(k^) is called an "optimal" mode, and its 
wavevector is called {Qlt}- 

The ideal case is that we can build a spin state satisfying 
two conditions 

Condition (1) {s^} are entirely linear combinations of 
optimal LT modes 

Condition (2) \si\ 2 = 1 everywhere (unit length con- 
straint) 

If both conditions are satisfied, these must be ground states, 
and all ground states must be of this form. 

In the case of a Bravais lattice (m = 1), the LT modes are 
just plane waves e ik r , and one can always construct a planar 
spiral configuration [1], s(r) = cos(Qlt ■ r)B + sin(QLT ■ 
r)C, where B and C are orthogonal unit vectors, and the spa- 
tial dependence consists only of optimal modes [2|. In the 
simplest cases, Qlt is at high symmetry points on Brillouin 
zone corners, and one can construct a combination of opti- 
mal modes which is ±1 on all sites, which defines a collinear 
ground state, as in the phase diagrams in Ref.[32l 

Thus, non-Bravais lattices are necessary in order to get non- 
coplanar states. (But not sufficient: it appears that, on non- 
Bravais lattices with high symmetry, the commonest ground 
states are still collinear or coplanar.) In lattices- with-a-basis, 
however, the LT eigenmodes have different amplitudes on dif- 
ferent sites within the unit cell, and it is not generally possi- 
ble to make any three-component linear combination of the 
best modes that satisfies the unit-length constraint. (There is 
an exception for lattices in which the neighbors -of-neighbors 
are all second neighbors, such as the diamond 11331 or honey- 
comb [34] lattices. 

A "generalized" L-T method for non-Bravais lattices was 
introduced by Ref.|35](see also Ref. 36) and applied to spinels 
with both A and B sites magnetic J2] [3T|. However, this 
method involves site-dependent variational parameters, so one 
must already understand the pattern of the ground state in or- 
der to make it into a finite problem; in practice, this method 
appears quite similar to our method (Sec.|V]and |A2[ ) of pro- 
jecting a layered structure to a one-dimensional chain. 

Although the LT optimal modes (usually) give the exact 
ground states in the cases we focus on, we believe the exact 
ground state is frequently built mainly from almost-optimal 
modes; that is, although a linear combination of optimum LT 
modes violates the unit-spin constraint, with a small distortion 
it may satisfy the constraint and be the ground state. (That 
distortion necessitates admixing other modes but with small 
amplitudes, since they carry a large energy penalty, according 



to (2.2 1.) In particular, we anticipate that (for incommensu- 
rate orderings) the true ordering wavevector lies in the same 
symmetry direction as the LT wavevector; and that the phase 
diagram for optimum LT wavevectors mostly has the same 
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topology as the actual phase diagram for ground states. Thus 
the LT modes can serve as a "map" for navigating the param- 
eter space of { Ji} and for understanding the ground state spin 
configurations. 

An important caveat is that almost all Q vectors have 
symmetry-related degeneracies, and the LT analysis is silent 
on how these modes are to be combined with different spin 
directions, so the specification of the actual spin configura- 
tion is incomplete. (An example is the "double-twist" state, 



of Sec. VI ) As a corollary, a single phase domain on the LT 



mode phase diagram might be subdivided into several phases 
in the spin-configuration phase diagram, that represent differ- 
ent ways of taking linear combinations of the same LT modes. 
This cannot be detected at the LT level. 

One immediate insight is afforded by considering the LT 
phase diagram. Shor t-range couplings have Fourier trans- 
forms J Q( g(k) in (2.1 1 that vary slowly in reciprocal space. 



Such functions typically possess extrema at high-symmetry 
points in the Brillouin zone; the same is probably true for the 
optimum eigenvalues and their wavevectors Qlt- That cor- 
responds to simple, commensurate ordering in real space. In 
order to get the optimal LT mode (and presumably the actual 
ordering) to be incommensurate, or to possibly stabilize states 
with stacking directions other than (100), one needs to include 
more distant neighbor couplings. 

In practice, we never used LT to directly discover the 
ground state spin configuration; its value is to quickly prove a 
given state is a ground state. But the LT viewpoint did inform 
the Fourier-transform diagnostic we used in analyzing the out- 



puts of iterative relaxation (Sec. II B Furthermore, when we 
operated in the "designer" mode (seeking the couplings that 
stabilize a specified state) we used the LT modes as a guide 
or clue: namely, we found the {Jj} that made the ordering 
wavevector Q of our target state to be the optimal Qlt, which 
is easier than making Q be the ordering wavevector of the ac- 
tual ground state. 



B. Iterative minimization 

Our prime tool for exploration was iterative minimization 
starting from a random initial condition. Random spins are 
selected in turn and adjusted (one at a time) so as to minimize 
the energy, by aligning with the local field of their neighbors, 
till the configuration converged on a local minimum of the 
Hamiltonian. [37]. (Our criterion was that the energy change 
in one sweep over the lattice was less than a chosen tolerance, 
typically 10~ 9 ). 

It might be worried that such an algorithm gets stuck in 
metastable states, unrepresentative of the ground state; such 
"glassy" behavior is indeed expected in the case of Ising (or 
otherwise discrete) spins, or in randomly frustrated systems 
such as spin glasses. However, vector spins typically have 
sufficient freedom to get close to the true ground state IT381 
[391 . The typical ways they deviated from the ground state are 
just long-wavelength wandering ("spin waves") or twists of 
the spin directions. 

The only problem with the dynamics is that our algorithm 



is a variant of "steepest descent", one of the slowest of relax- 
ation algorithms, such deviation modes are indeed slow re- 
laxing For this kind of (local) dynamics, the relaxation rate of 
a long-wavelength spin wave at wavevector q is proportional 
to |q| 2 , i.e. 1/L 2 for the slowest mode in a system of side 
L. (In future applications, some version of conjugate gradient 
should be applied to give a faster convergence, or - if there is 
a problem is finding the right valley of the energy function - 
one might adapt Elser's "difference map" approach to global 
optimization l40l .) 

For initial explorations, we usually used very small cubic 
simulation boxes of L 3 cells (L = 3, 4, or 5). For each set of 
{Jij} tested, we tried both periodic and antiperiodic bound- 
ary conditions, as well as even or odd L. Usually one of those 
four cases accomodates an approximation of the infinite sys- 
tem ground state, though of course any incommensurate state 
must adjust either by twisting to shift the ordering wavevector 
to the nearest allowed value, or else (as we observed) via the 
formation of wall defects. We tried to distinguish the ground 
states which were "genuine" in that a similar state would re- 
main stable in the thermodynamic limit. In particular, out of 
the four standard systems we tried (even/odd system size, pe- 
riodic/antiperiodic BC's), a "genuine" state should be the one 
with lowest energy. 

For a large portion of the parameter space, the ground states 
were planar spirals, essentially no different from the solutions 
guaranteed in the Bravais lattice case. Many of the noncopla- 
nar configurations found were "non-genuine" artifacts of finite 
size when the periodic (or antiperiodic) boundary conditions 
and dimension L were incompatible with the natural period- 
icity of the true ground state. One might expect the wrong 
boundary condition to simply impose a twist by ir/L per layer 
on the true ground state, but instead the observed distortion of 
the spin texture was sometimes of the natural periodicity of 
the ground state, a "buckling" occurred; that is, the configura- 
tion consists of finite domains similar to the true ground state, 
separated by soliton-like domain walls. 

The greatest difficulty in our procedure was not obtain- 
ing an approximation of the ground state, or even deciding 
whether it was genuine. Rather, it was grasping what the 
obtained pattern is, and how to idealize it to a periodic (or 
quasiperiodic) true ground state of the infinite system. We 
were aided by the following three diagnostics. 



1. Diagnostic: Fourier transform: 

Configurations obtained by iterative minimization were 
Fourier transformed and the norms of each Fourier component 
J2 S |s(k, s)| 2 were summed (combining the sublattices) fill . 
This suffices to identify the state when it is a relatively sim- 
ple antiferromagnetic pattern, or an incommensurate state de- 
scribed as a layer stacking. In any case, the results can be 
compared to the LT mode calculation to see if the found state 
achieves the LT bound. 
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FIG. 1: (COLOR ONLINE). Common-origin plots for four kinds 
of spin ground states on the octahedral lattice, appearing for dif- 
ferent choices of the exchange interactions {Ji}. (a), planar spiral 
(b). cuboctahedral state (c). asymmetric conic spiral (d). alternating 
conic spiral. 



three nonzero eigenvalues, respectively. (Rotating the spin 
configuration so the principal directions of M are the coor- 
dinate axes usually manifests the spin configuration's symme- 
tries.) 

C. Variational Optimization 

Through the diagnostic techniques described in the previ- 
ous subsections, it is normally straightforward to parameter- 
ize the spins in the ground as s(r, {a}), i.e. as a function 
of position and some arbitrary set of parameters, {a}. The 
exact values of these parameters can (formally) be calculated 
from optimizing H(s(r, {a}), {J}). One major advantage of 
relying upon this method is that it reduces the effect of numer- 
ical artifacts from iterative minimization. For example, we no 
longer enforce an arbitrary periodicity upon the LT wavevec- 
tor. 



2. Diagnostic: common-origin plot: 



D. Conceptual framework for bridging states 



The simplest visual diagnostic of a state is the "common- 
origin plot", in which each spin's orientation is represented as 
a point on the unit sphere. For example, an incommensurate 
coplanar spiral state would appear as a single great circle on 
the common-origin plot (see Fig. [TJ. 

A drawback of the common-origin plot is the lack of infor- 
mation on the spatial relation of the spins. (For example, a 
"cone" might appear consisting of closely spaced spin direc- 
tions, but these might belong to widely spaced sites.) Further- 
more, this diagnostic is quite fragile in configurations where a 
domain wall or other defect has been quenched in. 

In the case of one-dimensional chains (see Sec. |V B| i, we 
may instead use the "end-to-end" spin plot, where the tail of 
each spin vector is on the head of the previous spin vector. The 
advantage is that (i) images are not so obscured by overlaying 
of different vectors, and (ii) spatial information is captured, in 
particular defects where the state has "buckled". l42l 

In tandem with the common-origin plot, an algorithm was 
used that found groups of spins with (nearly) identical spin 
directions, within a chosen tolerance (typically, a dot prod- 
uct greater than 0.99). For commensurate patterns comprising 
a finite set of spin directions, this allowed the magnetic unit 
cell and spin configuration to be read off, but it was useless 
for incommensurate states, in which no spin direction exactly 
repeats. 



This subsection is not about a technique, but a classification 
of two ways that ground states may be related to each other, 
and thus of two kinds of phase boundaries in the phase dia- 



grams (Sec. VII and Appendix|A]i. We call these two concepts 
"encompassing states" and "families of degenerate states". 



1. Encompassed states 

We call a ground state "encompassed" if it is a special case 
of another, more general state. For example, a ferromagnet 
is encompassed by a helimagnet, since letting the helimag- 
netic angle go to zero produces a ferromagnet. "More gen- 
eral" means there is a continuous family of states such that 
each particular combination of couplings Jj's completely de- 
termines a particular member of that family. Moreover, the 
(more general) encompassing state necessarily spans at least 
one more dimension of spin space, so if the encompassed state 
is coplanar, the encompassing one is non-coplanar. However, 
while every encompassing state is more general than the state 
it encompasses, not every more general state will encompass 
a particular ground state. For example, the asymmetric conic 
is more general than the splayed ferromagnet (both of these 
states are defined in Sec. [V]), but it does not encompass the 
splayed ferromagnet; there is, however, a third class that en- 
compasses both these classes. 



3. Diagnostic: Spin moment-of-inertia tensor 
We computed the 3x3 tensor 



(notice Tr(M) = 1) and diagonalized it. We recognize copla- 
nar or non-coplanar spin states as those where M has two or 



2. Degenerate states 

By contrast, "degenerate states" means that for certain 
combinations of Ji 's, there is a continuous family of exactly 
degenerate ground states. Most commonly, this is the result 
of decoupling between sublattices of spins, meaning one can 
apply a global rotation limited to just one of the sublattices 
while remaining in the degenerate manifold. This can come 
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about in two ways. The trivial way is when all Jj's that cou- 
ple those sublattices vanish. The more interesting way is when 
the couplings are nonzero, but cancel generically in all the 
ground states; the simplest example of this kind is the J1-J2 
antiferromagnet on the square [3 1 or bcc [4| lattice, in the J2- 
dominated regime in which each of the even and odd sublat- 
tices realizes plain Neel order. Apart from decoupling, degen- 
erate manifolds are also sometimes realized by simultaneous 
rotations involving all sublattices with some mutual constraint 
(e.g. in the nearest-neighbor kagome lattice, the constraint 
that the spins add to zero in every triangle). 

Independent of the categories mentioned, these degenerate 
families may either be simply degenerate, meaning the ground 
state manifold is labeled by a finite number of parameters - 
one nontrivial angle in the case of the J1-J2 antiferromagnet 
— or else highly degenerate, meaning the number of param- 
eters scales with size as 0(L n ) with n > 0. An example with 
trivially decoupled, highly degenerate ground states is a lay- 
ered lattice with vanishing interlayer couplings, so the number 
of free parameters scales as O(L); the J = 1 kagome anti- 
ferromagnet is a mutally constrained, highly degenerate case 
with 0(L 2 ) parameters, i.e. extensively many. 

Both the simply and highly degenerate families also have 
clear signatures in reciprocal space. In the simple case, the op- 
timal wavevectors are discrete and related by symmetry, e.g. 
the J1-J2 square lattice antiferromagnet has Qlt = (1/2, 0) 
or (0, 1/2); in this case, the degeneracy lies in the freedom 
to mix these degenerate LT eigenmodes with different coef- 
ficients, not from the presence of extended (one-, two-, or 
even three-dimensional) surfaces in the Brillouin zone. By 
contrast, in the highly degenerate scenario, typically holding 
for special combinations of couplings, the LT optimum wave- 
fectors occur not just at isolated Qlt in the Brillouin zone, 
but on extended (one-, two-, or even three-dimensional) sur- 
faces, i.e. one has degenerate eigenmodes of the LT matrix 
that are not symmetry-equivalent. The rhombohedral lattice 
with J 1, J 2, and J 3 has a degenerate one-parameter fam- 
ily of wavevectors corresponding to different coplanar spi- 
rals 1143 1 . The three-dimensional pyrochlore lattice with only 
nearest-neighbor (Ji) couplings is a well-known example of 
the highly degenerate scenario, requiring extensive number 
of parameters. In that case, the minimum LT eigenvalue is 
uniform throughout the Brillouin zone 1441 (a so-called "flat 
band"). 



3. Encompassed and degenerate states as bridges in phase 
diagram 

What encompassing states and families of degenerate states 
have in common is to serve as bridges between simple states. 

In the "encompassing" case, the encompassing state is typ- 
ically stable in a domain of parameter space of nonzero mea- 
sure. When one adjacent phase in a phase diagram is encom- 
passed by the other, they are necessarily related (in our T = 
phase diagram) by a continuous transition, usually involving 
a symmetry breaking. 

In contrast, degenerate families are (frequently) confined 



to phase boundaries. Even when they occupy a finite area in 
a slice of parameter space (e.g. the ( Ji, J2) plane when all 
other couplings are zero), turning on additional couplings can 
remove the degeneracy. 

A corollary is that the naive classification of continuous 
or first-order phase transitions does not work. Consider two 
phases separated by a phase boundary on which a degenerate 
family is stable. Each of the two phases (or the limit of ei- 
ther as the boundary is approached) is a special case from the 
degenerate family. Since the limits taken from the opposite 
directions are different, it appears at first as an abrupt tran- 
sition. On the other hand, it is possible to take the system 
continuously from one phase to the other if we pause the pa- 
rameter variation when we hit the phase boundary, and follow 
a path through the degenerate manifold from one of the limit- 
ing states to the other one. 

Furthermore, turning on additional parameters generically 
destroys the degeneracy. That converts the degenerate family 
into an encompassing family, and the single phase boundary 
into two continuous ones. Specifically, starting in one of the 
main phases, we cross a small strip of phase diagram in which 
the configuration evolves (determined by the parameter com- 
bination) from one of the limiting states to the other one, and 
then enter the other of the main phases. Thus, the "encom- 
passed" kind of transition is distinct from either a first-order 
transition (between two unrelated states) or an ordinary con- 
tinous one, and will be indicated on phase diagrams with a 
distinct kind of line. 



E. Cluster analysis: two degenerate ground states 

The "cluster" method is a rigorous analytic approach to 
ground states, alternative to the LT mode approach. 11451 . 
which depends on decomposing the Hamiltonian into terms 
for (usually overlapping) clusters, and finding the ground 
states for one cluster. If these ground states can be patched to- 
gether so as to agree where they overlap, the resulting global 
state must be a ground state and all ground states must be de- 
composable in this fashion. In this way we can characterize 
the degenerate states appearing for two special combinations 
of Jj's. 



1. Antiferromagnetic Ji only 

In this case, the cluster is a triangle (one face of an octahe- 
dron, including one site each from the x, y, and z sublattices). 
The ground state of such a triangle is the usual 120° arrange- 
ment of spins. If all such triangles are to be satisfied, then 
wherever two of them share an edge, the respective unshared 
spins are forced to have identical directions - in the present 
case, spins on opposite corners of the octahedron. Thus, a line 
of x sublattice spins in the x direction (or similarly of the other 
sublattices in their directions) is constrained to be the same. 

This high degeneracy is not limited to the single point in 
parameter space J\ < 0. If we turn on J 2, which couples the 
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nearest neighbors aling those lines of spins, the same config- 
urations remain the ground state until J 2 is negative and its 
magnitude sufficiently large compared to | JjJ; less obviously, 
the same thing is true for J 3 = J4, varied together. 

This allows two different kinds of highly degenerate state: 
(a) One sublattice (say x) has = +A along every line. 
Within the other sublattices, each yz plane has an independent 
rotation about the A axis. Thus the spin directions are 



s 1 = -1a±^B(x), 



(2.4) 



where B is a different unit vector in each plane, and we take 
the + or — sign in the y or z sublattice, respectively. 

The common-origin plot for this state looks superficially 
like a conic spiral, the cone being formed by the y and z spin 
directions. In reality, whereas an incommensurate spiral gives 
a uniform weight along the spiral in the common-origin plot, 
this state gives a random distribution which approaches uni- 
formity only in the limit of a very large system. 

(b). For a second family of (discretely) degenerate states, 
we choose 

s( ni + 1/2, n 2 ,n 3 ) = i[fJ, / 2 (n 2 ), -/ 3 (n 3 )];(2.5a) 



s{n 1 ,n 2 + 1/2, n 3 ) 
s(n 1 ,n 2 ,n 3 + 1/2) 



V2 
1 

1 



hAK),0,/ 3 (n 3 )];(2.5b) 
[AK),-/ 2 (n 2 ),0].(2.5c) 



where fi(rii) = ±1 are arbitrary. Notice that (2.5c 1 uses (a 



subset of) the cuboctahedral directions. Typically, in a suf- 
ficiently large system, all those directions are used nearly 
equally; the common-origin plot would show a cuboctahe- 
dron. However, the spins do not have a regular pattern in space 
since (2.5c 1 is random, with a d iscret e degeneracy O(L) in a 
system of L 3 cells. The states (2.5c 1 represent a degenerate 



family of states, as formulated in Sec. II D the optimum LT 
eigenvalues are found at all wavevectors Q lying on the (100) 
axes. 

We are not interested in the high degeneracy for its own 
sake; its significance is that various kinds of ordered states 
can be selected out of it, by turning on additional couplings 
(even infinitesimally). Thus, the high-degeneracy parameter 
combinations will be corners of phase domains in the phase 
diagram. 



2. Antiferromagnetic J± and J2 



Let L Q be the net spin of the octahedron centered on a. 



(2.6) 



where a is a cubic lattice vertex and i <E a means site i is on 
one of the six bonds from a, Consider a Hamiltonian written 

as 



J 



On the one hand, expanding the square shows this is sim- 
ply the antiferromagnet with J\ = J 2 = J. On the other 



hand, it is obvious from (2.7i that any configuration with a 



net (classical) spin of zero on every octahedron is a ground 
state. This is another example of a degenerate ground state 
family (Sec. |IID) ; in this case the continuous degeneracy is 
macroscopic. This Hamiltonian is constructed in exactly the 
same way as those of well-known highly frustrated lattices 
(kagome, checkerboard, half-garnet, pyrochlore) that have 
similar ground state degeneracies. 



III. COPLANAR STATES 

Several different collinear or coplanar ground states can be 
stabilized within the octahedral lattice, only one of which re- 
quires couplings beyond J 2 . We will describe them from the 
smallest to the largest magnet unit cells. The most elemen- 
tary of these is the ferromagnetic state, in which all spins are 
aligned in the same direction, and which obviously requires 
predominantly positive couplings. This state is composed of 
a (0,0,0) LT mode with equal amplitudes on every sublattice 
(so the normalization condition is already satisfied). 

There is also the "three-sublattice 120° antiferromagnetic 
state", whose unit cell is the primitive cell. Each of the three 
sublattices has a uniform direction; the net inter-sublattice 
couplings are antiferromagnetic, so (as in the ground state 
of a single antiferromagnetic triangle) the respective spin di- 
rections are 120° apart and coplanar. Thus this state, too, is 
characterized by ordering wavevector Q = (0, 0, 0), but not 
the same LT mode as the ferromagnetic state. Instead, this 
one is from the two degenerate modes at Q = (0, 0, 0) that 
are orthogonal to the uniform mode. (Any combination of 
these modes has unequal magnitude on the different sublat- 
tices, which is why both modes need to be present in the spin 
state, combined with different spin directions.) This state is 
a special case of the highly degenerate ground states found 
when only J x < (Sec.[rTEl). 



(2.7) 



The next group of coplanar states are the Q / antiferro- 
magnetic states, of which there are three kinds, characterized 
by having ordering wavevectors of type (0,0,1/2), (0,1/2,1/2), 
or (1/2,1/2,1/2). Each of these states is antiferromagnetic 
overall within every sublattice; the sublattices decouple, since 
any inter-sublattice interaction couples a spin in one sublattice 
to equal numbers of spins pointing in opposite directions in 
the other sublattice. These are LT states; in the (0,0,1/2) and 
(0,1/2,1/2) cases, the LT mode used is nonzero on only one 
sublattice, and a different one of the three symmetry-related 
wavevectors is used for each sublattice (the one with the same 
distinguished direction), e.g. the x sublattice uses (1/2,0,0) 
or (0,1/2,1/2) modes. Notice that in these two cases, the 
spins repeat ferromagnetically along some directions (within 
a sublattice); this is a consequence of the anisotropy of the 
intra-sublattice couplings. Qualitatively, these states are sta- 
ble when J 2 is different from J 2 . All of these states can be 
realized with collinear spins. 

Lastly, the octahedral lattice admits helimagnetic states. 
These states require at least J 3 couplings to become stabilized. 
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The helimagnetic states must be composed entirely out of 
(q, q, q) modes. This is because any other wave-vector would 
break the symmetry between the sublattices. More precisely, 
helimagnetic states are generically a function of one variable, 
k ■ r, and are therefore equivalent to a one-dimensional sys- 
tem. It will therefore be amenable to stacking vector analysis, 
developed in [V] But using stacking vectors to transform the 
octahedral lattice to a one-dimensional chain will necessarily 
produce a non-Bravais lattice unless the stacking vector (111) 
(or a permutation of sign). And any helimagnetic mode in the 
one-dimensional non-Bravais lattice will break normalization 
in the octahedral lattice, since some spin directions would be 
represented more than others (this will be allowable for con- 
ies because they mix multiple modes, but helimagnets are ex- 
plictly single mode). Therefore, the only allowable stacking 
vector (and by implication, wave-vector) is (1 1 1). 



IV. CUBOCTAHEDRAL STATES 

The octahedral lattice possesses two kinds of "cuboctahe- 
dral" state, stable in different domains of parameter space, for 
which the common-origin plot takes the form of a cuboctahe- 
dron, i.e. twelve spin directions of the form (1, 1, 0)/V2 and 
its permutations [Fig.[TJb)]. The magnetic unit cell is 2 x 2 x 2 
for both of these true cuboctahedral states (spuriously cuboc- 
tahedral states were remarked in Sec. |IIE \ \ . They differ in 
that the angles between neighboring spins (which are in dif- 
ferent sublattices) is 60° in one kind of cuboctahedral state but 
is 120° in the other kind. 

As worked through in this section, the cuboctahedral states 
can be understood from any of three approaches: 

(a) Cluster construction: the Hamiltonian can be decom- 
posed into a sum of terms, each for an octahedron; we 
can patch together the ground states of the respective 
octahedra to obtain a ground state of the whole lattice. 
(For 60° cuboctahedral only.) 

(b) Degenerate perturbation theory: two special sets of cou- 
plings give degenerate families of ground states, out of 
which a small additional coupling can select the cuboc- 
tahedral state. (For 120° cuboctahedral only.) 



(c) The Luttinger-Tisza framework of Sec. |II A| (For both 
kinds of cuboctahedral state.) 

The 120° cuboctahedral state is a subset of the Jj-only 
antiferromagnetic (Ji < 0) ground states described above 
(Sec. |IIE 1 1 > - Thus, this state is stabilized even in the limit 
as J2 (or other distant couplings) become arbitrarily small. It 
is the only non-coplanar state we found that does not require 
any couplings beyond J 2 and J' 2 . 



A. Lattice as union of cuboctahedral cage clusters 

The first cuboctahedral state noticed was in the J\-J-2 mag- 
net on the kagome lattice [46 47 1, with J\ ferromagnetic and 



J2 antiferromagnetic. There is a range of ratios J2/J1 in 
which the magnetic unit cell on the kagome lattice is 2 x 2. 
Taking that cell as given, the possible ground states are those 
of the twelve-site cluster made by giving periodic boundary 
conditions to one unit cell - a cluster which is topologically 
equivalent to a single cuboctahedron (even when couplings to 
any distance are taken into account), and hence include the 
cuboctahedral state. 

Turning to the present case of the octahedral lattice, in fact 
this is a union of cuboctahedral cages surrounding cube cen- 
ters, complementary to the octahedral clusters surrounding the 
cube vertices. We can apply the "cluster" construction (see 



Sec. II E 1 to these cages by representing its lattice Hamiltonian 
as a sum of cuboctahedron Hamiltonians, with ji = Ji/2, 

3i = ^2/2, 32 = J'2/^^ an d ji = ^4/2 (Ji, J2, J'2 an d J4 are 
shared by two cuboctahedra); also j'3 = J3 and j\ = J{. 

Consider for a moment the ground state of 12 spins placed 
on an isolated cuboctahedron, [48 1 with couplings j±, fa, fa, 
and fa. Now, it is well known that, on a chain (i.e. a dis- 
cretized circle), fa > and (small) 32 < give a gradual, 
coplanar spin spiral; on a circle with the right number of sites 
the spin directions point radial to that circle. Roughly speak- 
ing, the three-dimensional analog of this happens on a cuboc- 
tahedron, which is a discretized sphere: if fa is ferromagnetic 
and one of the more distant couplings is antiferromagnetic, 
the spin ground state is a direct image of the center-to-vertex 
vector in the cuboctahedron (modulo a global 0(3) rotation of 
the spins). This notion only works for the 60° kind of cuboc- 
tahedral, occurring for J\ > 0, in which the nearest neighbors 
are (relatively) close to being parallel. 

To build a global state in which every cuboctahedron has 
the cuboctahedral spin configuration, the difficult part is just 
to make the spins agree where they are shared between cuboc- 
tahedra: that is achieved by applying mirror operations in al- 
ternate layers of cuboctahedra, such that e.g. the sf compo- 
nents are multiplied by (— l) x / a . 



B. As special case of Ji -only antiferromagnet 

When we have only J\ < couplings, the ground state is a 
degenerate state with 120° angles between nearest neighbors, 



as written in Eq. ( 2.5c 1 of Sec. II E 1 In that degenerate state 
of Sec. |II E 1 1 the spins take (some of or all of) the cubocta- 
hedral directions, but do not have a genuine cubic spin sym- 
metry. As soon as an arbitrarily small antiferromagnetic J 2 
is added as a perturbation, a subset of these states is selected, 
which is the 120° kind of cuboctahedral state. In the notation 
of ( |2.5c[ ), this true 120o cuboctahedral state takes the follow- 
ing form: 

s(»i + ~, n 2 ,n 3 ) = i[0, (-1)" 2 , (4.1a) 

s(m,n 2 + *,n 3 ) = i[(-ir,0,±(-l)^];(4.1b) 

s(n 1 ,n 2 ,n 3 + ~) = ±-^[(-ir,(-l)" 2 ,0].(4.1c) 
where + corresponds to the 120° state and — to the 60°. 
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C. Luttinger-Tisza approach to cuboctahedral states 

Alternatively, both cuboctahedral states can be understood 
within the LT framework. For certain domains of parame- 
ter values, the optimal LT modes have wavevectors Qlt of 
{1/2,0,0} type. It can be worked out that for e.g. Qlt = 
(1/2, 0, 0), one eigenmode has amplitudes (1, 0, 0) on sublat- 
tices (x, y, z), i.e. its support is only on the x sublattice. Each 
of the other two eigenmodes has its support equally on the y 
and z sublattices, the amplitudes being (0, 1,±1). When the 
first eigenmode is optimal, we get the decoupled (1/2, 0, 0) 



antiferromagnet already described in Sec. Ill when either of 



the two-sublattice eigenmodes is stable, a cuboctahedral states 
is found. 

It is obviously impossible to satisfy normalization in a spin 
state using just one of the two-sublattice modes, since its am- 
plitude vanishes on the third sublattice. To build a normal- 
ized ground state, it is necessary and sufficient to form a lin- 
ear combination using all three of the symmetry-related Qlt 
wavevectors, associating each with a different orthogonal spin 
component. Thus the spin directions are (1, 1, Q)/v2, with 
all possible permutations and sign changes, as we already 
saw in Eq. (4.1 1. These states could be called a commensu- 



rate triple-Q state. The eigenmode with amplitudes of form 
(0, 1, 1) gives the 60° cuboctahedral whereas the one of form 
(0, 1, —1) gives the 120° cuboctahedral state. 



1. Absence of (1/2, 1/2, 0) Cuboctahedral State 



decoupling of the J2-only simple cubic antiferromagnet.) In 
light of these decouplings, we cannot call this state a (1/2, 1/2, 
0) cuboctahedral; it is merely a particular configuration out of 
a degenerate family that also includes collinear states. 



V. CONIC SPIRAL STATES 

The conic spiral states are generically incommensurate and 
constitute the most common class of non-coplanar state that 
we found. They are layered states, where the spins are all 
parallel in a given layer. That is, the lattice breaks up into lay- 
ers, normal to some stacking direction Qstack in real space. 
We encountered only Q s tack = {Q00} stacked conic spirals, 
so we concentrate on that case, but stacking directions other 
than {100} should be feasible in principle. Within each layer, 
all the spin directions are the same; as you look in each suc- 
cessive layer, the spin directions rotate around an axis c in 
spin space. Due to this layering, it is possible to map a conic 
state to a one-dimensional "chain lattice" (as introduced in 
Sec. |VB) , which is a significant simplfication in the analysis. 

Considered from the LT viewpoint, a conic spiral is a mix of 
two different modes: two spin components follow an incom- 
mensurate wavevector Q, spiraling as in a helimagnet, and 
the third component follows a commensurate wavevector Q' 
(Q' = 0[AF] in the asymmetric -conic case). In both Q and 
Q', all components transverse to the stacking direction must 
be zero. 



From the LT viewpoint, one would naively expect to con- 
struct similar noncoplanar cuboctahedral states of cubic sym- 
metry using Qlt of (1/2,1/2,0): why are they absent? After 
all, the {1/2, 1/2, 0} type wavevectors are threefold degener- 
ate, just like the (1/2,0,0) wavevectors from which the cuboc- 
tahedral states are built, and it is straightforward to follow the 
analogy of those states to construct a (1/2, 1/2,0) cuboctahe- 
dral (just allotting each mode one of the three cartesian direc- 
tions in spin space). Furthermore, if we include J4 and J' 4 
couplings in the Hamiltonian, there is a certain region of pa- 
rameter space in which Qlt = (1/2, 1/2, 0) can indeed be opti- 
mal, with the optimal LT eigenmodes being orthogonal to the 
eigenmodes that make up the (1/2, 1/2, 0) type antiferromag- 
netic state. Hence in that region, the putative (1/2, 1/2,0), 
cuboctahedral state really is a ground state. 

But closer examination of the LT matrix for k = 
(1/2, 1/2,0) shows that neither the cuboctahedral state, nor 
any noncollinear state, is forced. At this high symmetry point 
in the Brillouin zone, all intersublattice contributions to the LT 



matrix J a p(k) [see Eq. (2.1 1] cancel. That means the (J a /3) is 



a diagonal matrix. Its eigenvalues are —J 2 + 2 J 2 — 4 J4 — 2 J 4 
for the mode of the (1/2,1/2,0) antiferromagnet, plus two de- 
generate eigenvalues J 2 + 2J' 4 for the modes of interest here. 

Furthermore, the fact that only J 2 and J' 4 enter the formula 
indicates that all other couplings cancel out. Not only are 
spins of different sublattices decoupled, but each sublattice 
decouples into two interpenetrating (and unfrustrated) tetrag- 
onal lattices . (The latter decoupling is reminiscent of the 



A. Categories of conic spiral states 

The conic spirals divide into subclasses, the alternating and 
asymmetric conic spirals, according to whether they are sym- 
metric under the (spin space) symmetry of reflecting in the 
plane normal to c. In the octahedral lattice, the only type of 
conic spiral we observed was the asymmetric conic. However, 
in the chain lattice we find that two distinct classes of conic 
spiral are possible: the asymmetric conic and the alternating 
conic. In principle, both states should be possible in the octa- 
hedral lattice, but an analysis in terms of stacking vectors (see 



Sec. |V B| > reveals that longer range couplings (J 5 , J 6 , . . . ) 
would be necessary to stabilize the conies. 



1. Asymmetric conic spiral 

In the asymmetric conic, the ground states are linear com- 
binations of helimagnetic and ferromagnetic modes. Let the 
stacking direction be along the z axis, so all spins are only 
functions of z. (Alternatively, we could interpret z as the po- 
sition in the one-dimensional chain lattice of Sec. |V B| below.) 
The helimagnetic part is parametrized with a rotating unit vec- 
tor: 



A(z) = cos(Qz)A + sin(Q2:)Bo 



(5.1) 
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where (Ao, Bo, C) form an orthonormal triad. Then 

s(z) = cos(a)A(z) + sin(a)C; 



(5.2a) 



cos/3A(z + -) -sin(/3)C, (5.2b) 



with < a, f3 < ir/2. The spins in each sublattice rotate 
about some common cone axis C in spin space; spins of the 
same sublattice have the same component along the cone axis 
(giving net magnetic moments for both sublattices of N *sjj , 

(a) 

where is the component of a spin of sublattice a along 
the common axis). The different sublattices are antiferromag- 
netically coupled, so these net moments have different signs. 
And because the couplings within the sublattices are not equal 
to those of the other sublattice, the magnitude of the net mo- 
ments are not equal. This can be easier to understand if we 
think in terms of common-origin plots (see Figure ]T} , since 
then the spins all lie upon the surface of a sphere. In the com- 
mon origin plot, each sublattice forms a cone. These cones 
are along the same axis, but oppositely oriented and their az- 
imuthal (conic) angles are not equal. 



2. Alternating conic spiral 

In the alternating conic, one sublattice is a planar spiral, 
while the other is always a combination of the same helimag- 
netic mode and an antiferromagnetic mode. This state is thus 



represented, again using (5.1 



s(z) = cos(a)A(z) + (-l) z sin(a)C; (5.3a) 



A(z+2-); 



(5.3b) 



where < a < ir/2 for the alternating conic. Returning 
again to our common-origin plot (see Figure]!}, one sublattice 
forms a great circle along the equator of the sphere. The other 
sublattice now traces cones on each side of this equator (the 
common axis is the vector normal to the circle). The spins of 
the second sublattice alternate between the two sides of the 
equator, giving the antiferromagnetic component. Thus, the 
difference between sublattices is more fundamental in the case 
the alternating conies than in asymmetric conies. 



3. Splayed States 

There are two other states essentially related to these conic 
spirals, but are important enough to deserve their own names 
(this is much the same as ferromagnetism being a special case 
of helimagnetism). We term these states ferromagnetic and 
ferrimagnetic splayed states. 

Consider the alternating conic in the limit of the polar an- 
gle going to or 1/2. In both cases, the spins are confined 
to a plane, but they are emphatically not in a helimagnetic 
configuration. The sublattice that was helimagnetic is now 
ferromagnetic and the sublattice that was conic now merely 



alternates (that is, reflects about equatorial plane without ro- 
tation). If the polar angle is (1/2), then the dot product of 
spins in different sublattice is positive (negative) and the state 
is a ferromagnetic (ferrimagnetic) splayed state. 

While the difference between ferromagnetic and ferrimag- 
netic splayed state seems rather trivial here, it is more dra- 
matic when we think of asymmetric conies. The ferrimagnetic 
splayed state is produced when one of the conic angles and 
the polar angle go to 1/2 while the other conic angle remains 
arbitrary. But because the coupling between sublattices is an- 
tiferromagnetic for the asymmetric conic, the conic angles of 
the sublattices will always confine spins to opposite sides of 
the "equator." This means that the asymmetric conic will never 
continuously transform into a ferromagnetic splayed state, and 
so such a transition would necessarily be first order. 



B. Stackings and chain mapping 

For the cases of incommensurate conic spirals, our main an- 
alytic method is variational: we assume a functional form for 
the spin configuration (based on interative minimization re- 
sults) depending on several parameters, and optimize exactly 
with respect to them. Say we know that the correct ground 
state is stacked a stack of planes with identical spins — in 
practice this is determined empirically from the outcome of 
iterative minimization (Sec. II B i - then the variational prob- 



lem is equivalent to a one-dimensional (and hence simpler) 
one: layers of the 3D lattice may be mapped into a chain con- 
taining inequivalent sublattices. 

This mapping is fruitful in two or three ways. First, we 
could (and did) explore the chain lattice ground states using 
interative minimization, in much longer system lengths than 
would be practical in an L x L x L system. Second, it is uni- 
fying, in that various stackings of various three-dimensional 
lattices map to the same chain lattice. Finally, it illuminates 
what conditions are necessary in order to obtain non-coplanar 
states. 

Notice that if a stacked state is the true ground state of 
the three-dimensional lattice, its projection must be the true 
ground state of the chain projection (since the chain lattice 
states correspond exactly to a subset of three-dimensional 
states); but of course the converse is false: the proven opti- 
mal state of the chain lattice might be irrelevant to the three- 
dimensional lattice, when a different (e.g. unstacked) kind of 
ground state develops a lower energy. As coupling parame- 
ters are varied, that different state might become stable in a 
first-order transition; our only systematic ways to address that 
possibility are (i) iterative minimization (ii) watching for an 
exchange of stability between two LT eigenmodes at different 
wavevectors. And even with this method, we had to rely pri- 
marily upon iterative minimization for reliable results, as LT 
analysis is insufficient to determine the ground state, particu- 
larly in cases where a ground state cannot be constructed from 
the optimal LT modes. 

For the cases that concern us here, the chain lattice has a 
basis of two sites per cell, with inversion symmetry at each 
site; we take the lattice constant to be unity. The mappings to 
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chain sites z is given by 



Octahedral 



Q 



stack ' F 



(5.4) 



where Qstack is a vector of integers, having no common fac- 
tor. We let "even" sites be those with z integer and "odd" sites 
be those with z — integer+1/2. As in three dimensions, we 
consider inter-sublattice couplings j\ and j 3 , as well as intra- 
sublattice couplings out to distances 1 and 2, namely j 2 and 
j4 (between even spins) or j 2 and j' A (between odd spins). No- 
tice that, if j 2 = j'2 an( J J4 = i4> the chain system reduces to a 
Bravais lattice (with lattice constant 1/2) and its ground states 



are (at most) coplanar spirals, as explained in Sec. II A , that 
rules out Qstack = (1,1,1)- Notice that for stackings in low 
symmetry directions (and thus requiring larger coefficients in 
the Qstack vector), short range couplings in the octahedral 
lattice map to long range j'y couplings in the chain lattice, 
e.g. Qstack = (211) maps J x through J 4 to j x through j e . 
Because longer range couplings quickly appear, it is reason- 
able to explore them in the chain lattice. In order to organize 
our exploration of parameter space, we shall call ji, j 2 , and 
j 2 "primary" couplings; j'3, j'4, and j 4 are "secondary" cou- 
plings, and if necessary are assumed small compared to the 
primary couplings. 

We encountered Qstack = (1, 0, 0) stackings often enough 
in the iterative minimization, and we searched for Qstack = 
(1,1,0) type stackings also (although this search was ulti- 
mately unsuccessful). The Qstack = (001) mapping is il- 
lustrated in real space in Fig. [2] The numerical values of the 
mapped couplings is given by a matrix multiplication: 
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(5.5) 



where jo is the energy within a plane of constant z in the octa- 
hedral lattice. Notice that, for this stacking vector, couplings 
J 1 through J4 are projected down to jo through j 2 . This helps 
explain the absence of stable conies in the octahedral lattice. 
We only find conies in the chain lattice for j'3 or longer-range 
couplings. But for this stacking Qstack = (001), require at 
least a J5 coupling to generate an analogous coupling (for the 
alternating conic, J§ is more likely, given the asymmetry be- 
tween the sublattices). 



C. Transversely modulated spirals 

This is a hypothetical (but likely) class of states. Here 
"transversely modulated" means that when we decompose the 
lattice as a stacking of layers, a single layer does not have a 
single spin direction, but instead a pattern of spin directions. 
Whereas the asymmetric conic spiral used ordering wavevec- 
tors (say) (Q, 0, 0) and (0, 0, 0), and the alternating conic spi- 




FIG. 2: (COLOR ONLINE) The mapping between the octahedral 
and chain lattices for a (001) stacking direction. All sites with the 
same value of z on the x and y sublattices (red x's) are projected 
onto the same point, while those in the z sublattice (blue squares) are 
mapped onto a different point (again a function z). As such, there are 
twice as sites mapped onto any x as there are mapped onto a square. 
The spatial structure is shown in dotted lines, while couplings are 
shown in solid lines. 



ral used (Q,0, 0) and (1/2,0,0), a transversely modulated 
spiral might replace the first wavevector by e.g. (Q, 1/2, 1/2). 

Equivalently, if we look at a column of successive cells 
along the stacking direction, in a plain conic spiral (whether 
alternating or asymmetric), adjacent columns are in phase, but 
in the transversely modulated conic spirals, different columns 
are offset in phase according to a regular pattern. It should be 
possible to generalize the chain mapping to such cases, but we 
have not tried it. 



VI. DOUBLE-TWIST STATE 

Here we describe the incompletely understood "double- 
twist" state, which has attributes in common with both cuboc- 
tahedral and conic states, and was observed only for a small 
set of couplings (Ji = —2, J 2 < 0, J3=l, all others zero). 
These couplings were selected to give Qlt = (Q, Q,0), as 
the ground states previously encountered had a Qlt either 
along the (111) direction or on the edges of the Brillouin zone. 
For a given Q, the value of J 2 is determined by: 

J 2 = 4%/2 • (3cos 2 (Q/2) - l)/y/l - 2cos 2 (Q/2). (6.1) 

We particularly studied the Q = 3/8 case. (Note that iterative 
minimization necessarily probes commensurate states, due 
to our boundary conditions.) This corresponds to J2/J3 = 
—3.7717, according to Eq. ( |6.1| l. 

The double twist state is, to good approximation, com- 
posed solely of (Q, Q, 0) modes (for normalization, there 
will necessarily be other wave-vectors, but these have rela- 
tively small amplitudes). Unlike previous states, each sublat- 
tice has nonzero contributions from all (Q, Q, 0) wavevectors, 
rather than a subset. The weight of each sublattice in a given 
(Q, Q, 0) mode differs between sublattices, approximately in 
proportion to relative weight in the LT optimal mode with a 
similar Qlt- 
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The spatial variation produced by this combination of 
modes is complicated. There is a stacking axis in real space, 
which we take to be z without loss of generality. Spin space is 
characterized by three orthonormal basis vectors: C defines a 
conic axis, around which the other two basis vectors A and B 
rotate as a function of z: 

A(z) = cos(Qz)Ao — sin((5z)Bo; (6.2a) 
B(z) = sin(Qz)A + cos(Qz)B ; (6.2b) 
C = AO) x B(z) = A x B . (6.2c) 

We can parameterize this cartoon of the double twist state 



sin($ y )A(z) 
sin($ x )B(2) 
rcos($ a: ) cos($ y )Cy 
(6.3) 

with $x(c) = Qx — (f> x , where <fi x is an arbitrary phase, sim- 
ilarly & y (y) = Qy — 4> y - Note the coordinates r in Sj(r) are 
the actual sites for sublattice i, which are half-odd-integers in 
the i component. Thus, addition to the twisting of the basis 



VII. PHASE DIAGRAMS 
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vectors along the stacking direction, in (6.2i, there are spatial 
modulations transverse to the stac king direction that appear in 
the coefficients of A, B and C in (6.3 1 



While the form described by (6.3 1 is close to what we 
observe with iterative minimization, it unfortunately does 
not satisfy normalization: the ground state necessarily con- 
tains admixtures of non-optimal modes. (To satisfy normal- 
ization using only the (Q, Q, 0) modes would require four- 
component spins.) In Eq. ( |6.3| , (—a, b, b) should ideally be 
the amplitudes (on the three respective sublattices) of the LT 
eigenvector at (0, Q, Q), while T is a weighting factor that 



reduces the deviations of the spins in (6.3 1 from uniform nor- 
malization. 

What if we demanded, not normalization of all spins, but 
only that the mean-squared value of |S;(r)| 2 be one in each 
sublattice? Since each cosine or sine factor has mean square 
of 1/2, and since a 2 + 2b 2 = 1, we ought then to have To — 
r = ^/2. Projecting the actual result of iterative minimization 
onto such modes gave V w 1.36. Also, whereas b/a = 1.64 in 
the actual LT eigenvectors, we found b/a ~ 1.5 in the results 
of iterative minimization. 

The double twist state can be viewed as related to the hy- 
pothesized transversely-modulated conies or to the cubocta- 
hedral states. In particular, the composition of this state in 
terms of LT modes is more similar to the cuboctahedral states 
than to any other configuration. 

The LT mode underlying this state, according to \6.\\ , has 
a continuously variable wavevector as J 2 is varied. Due to the 
limitations of iterative minimization with periodic boundary 
conditions, we have not followed the evolution of the double- 
twist state; in particular, we do not know if it becomes incom- 
mensurate in both the stacking (z) direction and the transverse 
directions. 



To understand how the ground states outlined in Sec. |ffl|iV] 
|V] and VI fit together, it is necessary to examine the phase di- 
agrams of the octahedral and chain lattices. A series of rep- 
resentative cuts through the phase diagrams for both lattices 
gives us a general sense of their topology, and specifically in 
what regions non-coplanar states are stabilized. 

Of course, rescaling the couplings by any positive factor 
gives an identical ground state (with energy rescaled by the 
same factor). Therefore we present the phase diagrams in 
rescaled coordinates, normally Ji — > Ji/\J\\ (except when 
Ji=0). 

An important aspect of all the phase diagrams is the classi- 
fication of the transitions into first-order (discontinuous), en- 
compassed (continuous), or degenerate: the distinction be- 
tween the last two kinds was explained in Sec. |IID| When- 
ever a continuous manifold of degenerate states is found (al- 
ways on a phase boundary), it is labeled in the diagrams by 
"0(L d y representing how the number of parameters (needed 
to label the states) scales with system size. 



A. Octahedral Lattice 

In the octahedral lattice we are fortunate in that most kinds 
of states have energies that can be written exactly as a lin- 
ear combination of couplings (given in Table [I] 149]); the 
phase boundaries between such phases are simply the lines 
(more exactly hyperplanes) where the two energy functions 
are equal. Most other phase boundaries are handled analyti- 
cally, e.g. the helimagnetic state and its "encompassed" ferro- 
magnetic and antiferromagnetic states. The only phase bound- 
ary not determined analytically from a variational form was 
the double-twist state, for which we do not have an exact vari- 
ational form; in this case the boundary was approximation by 
the LT phase boundary. We would expect that approxima- 
tion to be accurate for any such complex phase that is built 
entirely from a star of symmetry-equivalent modes, provided 
the neighboring phase is built from other modes. 

Using this information, we can easily find the phase bound- 
aries of various states. To aid in graphical display, we 
will normalize all couplings by | J\ | and restrict attention to 
I ^3 1/| J\ I < 1- Phase diagrams will be plotted in the variables 
( J2, J 2) representing a slice with ( J\, J3) fixed. In all such 
slices, the second, third, and fourth quadrants of the phase 
diagram are dominated by antiferromagnetic phases of order- 
ing vector (1/2,1/2,0), (1/2,1/2,1/2), and (1/2,0,0) re- 
spectively. Recall that all of these are nontrivially decoupled 
states, in which distinct sublattices can be independently ro- 
tated due to cancellations of the inter-sublattice interactions. 
When thermal or quantum fluctuations are added to the de- 
scription, "order-by-disorder effects" J5] 51 typically select 
specific states from these manifolds that are collinear. The 
first quadrant is dominated either by the ferromagnetic phase, 
or (if Ji < 0) by the 120° 3 sublattice phase. Cuboctahedral 
phases may be found near the J' 2 — axis when J 2 > 0. 

The phase transitions are always first order in the octahedral 
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State 


Q 


Energy/spin 


— Ferromagnetic 


(000) 


-4Jx - Ji - 2J' 2 - 4J 4 - 24, 
Jx = Ji + 2 Ji 


3 sublattice- 120° 


(000) 


2Jx -Ji- 2 Ja - 4J 4 - 2 Ji 


(1/2, 1/2,0)- AFM 


(U,o) 


- J 2 + 2J 2 - 4J 4 - 2J 4 


(1/2, 1/2, 1/2)- AFM 


(- - i) 

I 2 ' 2 ' 2/ 


J 2 + 2 - 4J 4 - 2 Ji 


(l/2,0,0)-AFM 


(1,0,0) 


J 2 - 2 J' 2 - 4J 4 - 2 Ji 


7i73-Cuboctahedral 


(1,0,0) 


-2Ji - J 2 + 4J 3 + 4J 4 


27r/3-Cuboctahedral 


(1,0,0) 


2 Ji - J 2 - 4 J 3 + 4 J 4 


Helimagnet 


(qqq) 


-2Ji - 2J 4 - Ji - J|/(8Ji), 
J s = 2 Ji + J 2 + 2J 2 + 4J 3 
and Jl = 2 J 3 + 2 J 4 + J 4 



Ji = i,J 3 = o 



TABLE I: Ground States of the Octahedral Lattice, with ordering 
wavevector Q. 



lattice, with the following exceptions, which can be classified 
according to the three scenarios for bridging states outlined in 
Sec. |IID| (encompassing and degenerate). 

(1) The transition from the helimagnet to either the 
(1/2,1/2,1/2) antiferromagnet or to the ferromagnet is 
continuous, since the optimal wavevector varies con- 
tinuously along (q, q, q) until it hits the commensurate 
value ((1/2,1/2,1/2) or (0,0,0)), then stops; this is an ex- 
ample of an encompassing state. 

(2) Transitions between two antiferromagnet phases are al- 
ways degenerate, since the phase boundaries in param- 
eter space are given by J2=0 or J 2 =Q, which (triv- 
ially) decouple sublattices (implying a degenerate fam- 
ily of states). In the Brillouin zone, the wavevector can 
evolve continuously along (l/2,q,k) and (1/2, 1/2, q) (for 
the (1/2,0,0) and (1/2,1/2,0) antiferromagnets, respec- 
tively). 

(3) Lastly, transitions between states of the same Qlt are 
degenerate, occurring where two eigenvalues of the LT 
matrix for Qlt cross, as a function of changing param- 
eters. This is found for the (0,0,0) modes (ferromag- 
netism and the 120° state) and the (1/2,0,0) modes (both 
types of cuboctahedral states with each other and with 
the Qlt = (1/2,0,0) antiferromagnetic state). Be- 
cause the degeneracy is limited to different eigenmodes 
of the same wavevector, these are O(L ) degenerate 
transitions. 

Consider first the phase diagram produced with ferromag- 
netic Ji and no couplings beyond J 2l J' 2 (Figure [5J. In this 
case, we find only four states, all of them coplanar (these are 
outlined in|Hl|). What's particularly important here, though, is 
the way that the phase diagram divides up into four quadrants. 
This is a fairly generic feature that we will see in other phase 
diagrams. Along the J 2 = line between the (1/2,1/2,1/2) 
AFM and the (1/2,0,0) AFM, in all the phase diagrams, we 
get a degenerate (0(L 2 )) decoupled state, in which each J 2 - 
coupled line has an independent staggered spin direction. 




Ferro 



AFM 
1/2, 1/2, 1/2) 




L 5 -4 -3 



-1 1 2 3 

4/Ji 



FIG. 3: Octahedral Phase Diagram for Ji = 1, J3 = 0. Squares 
indicate couplings tested with iterative minimization. Solid lines 
denote first order transitions , dashed lines denote degenerate transi- 
tions (Sec. |IID) , where the scaling of degrees of freedom is labeled. 
Dotted lines indicate a second order transition from encompassing 
states. Regions shaded gray indicate a non-coplanar phase. 



Ji = i,j 3 = -i/io 
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1/2, 1/2, 0) 

O(L) 



Cuboc 
tt/3 

1/2,0,0) 
Ferro 



AFM 
1/2, 1/2, 1/2) |S 



Heli 

q,q,q) 

AFM 
1/2, 0, 0) 



L 5 -4 -3 -2 -1 1 2 3 4 5 

4/\Ji\ 

FIG. 4: Octahedral Phase Diagram for Ji=l, J 3 = -1/10. First 
instance of a non-coplanar state (the 7r/3 cuboctahedral state, in the 
shaded region) and of helimagnetism. 
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-1 



12 3 4 

4/\Ji\ 



FIG. 5: Octahedral Phase Diagram for Ji = — 1, J3=0. First in- 
stance of the 2n/3 cuboctahedral state. The large circle point is the 
case of only J\ couplings, corresponding to the (highly degenerate) 
120° state. 



Let's now examine how the ferromagnetic J\ phase dia- 
gram is modified by an antiferromagnetic J 3 (Figure ffl. First 
of all, it stabilizes the tt/3 cuboctahedral state, our first exam- 
ple of a non-coplanar phase. J 3 also stabilizes a (q, q, q) heli- 
magnet at the center of the phase diagram. The boundaries of 
this phase are quite sensitive to J 3 : as J 3 becomes more an- 
tiferromagnetic, the helimagnet's phase boundaries with fer- 
romagnetism and (1/2,1/2,1/2) antiferromagnetism move out- 
ward in opposite directions, so as to increase the region of 
parameter space that is helimagnetic. Meanwhile, the phase 
boundaries of helimagnetism with the (1/2,0,0) antiferromag- 
net and tt/3 cuboctahedral state move inwards in opposite di- 
rections, so as to decrease the region of parameter space that 
is helimagnetic. The result is that, as J 3 becomes more an- 
tiferromagnetic, the helimagnetic region of parameter space 
first grows and later shrinks until J 3 =— Ji/2 < 0, where it 
disappears entirely. 

Now we turn to the phase diagrams with antiferromagnetic 
J 1, first considering arbitrary J 2 , J 2 with J 3 = (Figure 
B). In this case, we still see the quadrant structure, at least 
qualitatively; the upper right quadrant now represents the (or- 
dered) three-sublattice 120° state. However, a strip between 
the upper quadrants is occupied by the 27r/3 cuboctahedral 
state, a non-coplanar state which only requires two non-zero 
couplings. Furthermore, we find the (highly degenerate) "Ji- 
only" state along the boundary between the 27i73 cuboctahe- 
dral state and the three-sublattice 120° state. It is interesting 
to note that at this point in parameter space, the type of tran- 
sition changes. Because of the vanishing of the couplings, the 
minimum eigenvalue is a constant along (q,0,0). This makes 
the phase transition at this point a degenerate ((^(L 1 )), con- 
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FIG. 6: Octahedral Phase Diagram for Ji = — 1, Js=l/2. The 
double-twist state is found in the small region near the origin. The 
striped region indicates this slice (J3 = 1/2) contains a phase 
boundary between the 120° 3-sublattice and ferromagnetic states; 
for the parameter set plotted, that region has states of the degener- 
ate (marked "O(L )" parameters) and nontrivially decoupled kind. 
If we varied J3 to pass through that region of parameter space, we 
would cross a degenerate phase transition. 



tinuous transition, even though it is elsewhere first order. 

For antiferromagnetic adding ferromagnetic J 3 (Figure 
|6j, we see the 27r/3 cuboctahedral state expand, changing the 
topology of the phase diagram (it now shares a boundary with 
the (1/2,0,0) antiferromagnet). In addition, near the triple 
point of the (1/2, 1/2, 1/2) antiferromagnet, (1/2, 0, 0) anti- 
ferromagnet and 27r/3 cuboctahedral state, we find the double 
twist state ll50ll . Note that the slice of parameter space shown 
here, J\ = — 2 J 3 , includes the phase boundary between the 
three-sublattice 120° state and the ferromagnetic state, and on 
the boundary has an extra degeneracy of the kind described in 
SecED] 

The double-twist can be understood as a selection from the 
family of degenerate states. Along the ,]' 2 — line (boundary), 
we have a degenerate form of decoupled state: sublattices de- 
couple by row, giving a large degeneracy of ground states. 
However, this decoupling depended upon the alternating or- 
der within each sublattice which led to cancellations. A spiral 
distortion (selected by adding J\ , J 3 to J 2 only state, with J 2 
sufficiently small) within a sublattice allows a non-cancelling 
interaction between it and another sublattice, which lowers the 
energy. 

The Ji < phase diagrams superficially resemble the the 
Ji > phase diagrams, with two different cuboctahedral 
states appearing around the J 2 > axis, and a helimagnet 
or double-twist state (respectively) appearing in a small edge 
below the phase diagram's center. 
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We have also considered the case of J\ = with an an- 
tiferromagnetic J3 (phase diagram not shown). This phase 
diagram, apart from the trivial change of normalizing the cou- 
plings by I J3I instead of | J\\, strongly resembles the case of 
antiferromagnetic J\ and J3 = shown in Figure [5]); the sole 
difference is that we now find the ir/3 cuboctahedral state in 
place of the 2n/3 cuboctahedral. 

Iterative minimization found asymmetric conic states along 
the phase boundary between the (1/2,1/2,1/2) antiferromag- 
net and the (1/2,0,0) antiferromagnet, but we believe these are 
artifacts, in the sense we will describe. This boundary cor- 
responds to LT modes degenerate over a plane of wavevec- 
tors, leading to a degenerate family of spin ground states with 
an arbitary wavevector. These are generically non-coplanar 
spirals, except the limiting states of this family are collinear 
"encompassed states" (in the nomenclature from Sec. II Di. 
Thus, although these conic spirals are valid ground states, we 
do not count them as non-coplanar, since that is not forced by 
the couplings. This is an instance where the overlap between 
"encompassing" and "degenerate" states is especially stark, 
as the family of degenerate states coincides with the class of 
encompassing states. 

The impossibility of forcing any conic spiral in the octa- 
hedral lattice is understood by using the mapping (5.5 1 of 



couplings from the octahedral lattice to the chain lattice (see 



Sec. |VB| i. We will see shortly (Sec. |VH B) that stabilizing ei- 
ther kind of conic requires a coupling j'3 or 34 in the chain 



lattice; for a (100) stacking vector, Eq. (5.5 1 takes octahedral 
coupings Ji through J4 to chain-lattice couplings jo through 
j'2, so clearly couplings J 5 , J 6 or longer are required (and suf- 
ficient) to truly stabilize conic spirals in the octahedral lattice. 



B. Chain Lattice 

The chain lattice ground states are significantly more com- 
plicated than those of the octahedral lattice. Analytically de- 
termining the optimal energy of even the helimagnet becomes 
difficult when couplings beyond j'3 are included. Therefore, 
while we can easily determine a variational form for the ener- 
gies, we cannot analytically determine the ground state when 
couplings j'3 or higher are introduced. Energies are given in 
table|Il] which are then numerically optimized to give the sub- 
sequent phase diagrams. We once again normalize by \ji\, but 
we now plot j/,/^! x j 2 /\ji\, rather than J 2 /\Ji\ X J' 2 /\-h\- 
This change of convention does not have great physical im- 
plication, as the difference between J 2 and J 2 in the octahe- 
dral lattice is distinct from the difference between j 2 and j' 2 . 
Lastly, by the definitions of the chain lattice, several proper- 
ties of the phase diagram follow immediately. First, simul- 
taneous exchange of j 2 with j' 2 and 34 with j' 4 will merely 
change the labeling convention to distinguish the two sublat- 
tices. The ground state in the chain lattice must therefore be 
invariant under this operation. Furthermore, when j 2 = j' 2 
and 34 = j' 4 , this exchange will not change anything. In this 
region of parameter space, there is no difference between the 
two sublattices and the chain lattice becomes a Bravais lattice 
with unit cell 1/2. From this fact and Sec. Ill A| it follows that 



State 


Q 


Energy 


Ferro 





-2ji - 32 - 32 ~ 2j 3 - J4 - j'4 


AFM 





2ji - 32 - j'2 + 2j:i - ji - j' 4 


Heli- 


2ip 


— 2ji cos ij) — (j'2 + j'2) cos 2%j) 


magnet 




—2js cos 3ip — (34 + ji) cos Aip 


splayed 


Oand 


(ji + Ja) 2 /2j2 + 31 - j'2 ~ ji - 34 


ferro 


1/2 




splayed 


Oand 


(jl+j-i) 2 /2j2+j2-j'2-ji-j'4 


ferri 


1/2 




Alter- 


1/2 


—2(ji cos ip + ]3 cos 3ip) cos a 


nating 


and 


+J2 — 34 — j'2 cos 2ip — j'4 cos Aip 


Conic 


2ip 


— (j2 [cos 2tp + 1] + 34 [cos Aip — 1] ) cos 2 Q 


Asym- 





— 2(ji cos ip + J3 cos Zip) cos a cos /3 


metric 


and 


+2(ii + 33) sin a sin pJ - j 2 - j'2 - ji - j'4 


conic 


2ip 


— (32 cos 2 a + j' 2 cos 2 /3)(cos 2ip — 1) 

— (ji cos 2 a + j' 4 cos 2 P)(cos Atp — 1) 



TABLE II: Parametrizations and Energies of the Chain Lattice. Q 
denotes the ordering wavevector(s), as a multiple of 2n. The en- 
ergy per unit cell is given. The splaying angle is given by cos a = 
-O'l + J3)/2J2 for the splayed ferromagnet or +(jj + j 3 )/2j 2 for 
the splayed ferrimagnet. 



states in this region are necessarily coplanar. 

To classify phase boundaries in the chain lattice, it is impor- 
tant to consider limiting cases (i.e. encompassed states in the 
nomenclature of Sec. |HD| i. States with variational parameters 
(the helimagnetic wave-vector q, as well as the conic/splay an- 
gles) will undergo a second-order transition when their param- 
eters reach a limiting value (0 or 1/2 for the helimagnet angle, 
or 1/4 for conic angles). Thus there is a second-order tran- 
sition between helimagnetism and either antiferromagnetism 
or ferromagnetism, as well as between alternating conies and 
every other state (except asymmetric conic). Asymmetric con- 
ies, on the other hand, have second-order transitions to ferro- 
magnetism, antiferromagnetism, helimagnetism, or ferrimag- 
netic splayed states, but not to ferromagnetic splayed states or 
alternating conies. The splayed states, meanwhile, can only 
have second-order transitions to ferromagnetism or antiferro- 
magnetism (depending upon which type of splayed state it is), 
to the other splayed state, or the appropriate conies. All other 
transitions are necessarily first-order. 

Consider first the case of ji ferromagnetic with no cou- 
plings beyond j 2 , j' 2 (Figure [7]). The phase diagram displays 
the same quadrant structure that we found in the octahedral 
lattice. However, the quadrant structure is not identical in the 
two lattices. First of all, the ground states are different in the 
chain lattice (helimagnetism and splayed states instead of var- 
ious forms of antiferromagnetism). Secondly, the topology of 
the first and second order transitions are reversed for the two 
lattices. Both of these phenomena can be explained by appeal- 
ing to the additional degrees of freedom in the octahedral lat- 
tice. Because the octahedral lattice has three spatial variables, 
it has ground states that cannot exist in the chain lattice. This 
includes families of degenerate states at the phase boundaries 
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FIG. 7: Chain Lattice Phase Diagram for ji = l, J3=0, J4=0, ^4=0. 
Note similarity of phase boundaries to Figure [3] despite the lack of 
other similarity. 



of the octahedral lattice, producing second order transitions 
(when there are second order transitions in the chain lattice, 
they are principally due to encompassing states). 

Next, we examine the case j% antiferromagnetic and ja fer- 
romagnetic (Figure [8]). Several states in the quandrant struc- 
ture are different from the ferromagnetic j-y case (i.e. fer- 
rimagnetic vs ferromagnetic splayed), but more interesting 
is the presence of the asymmetric conic - the first instance 
of a nonplanar state in the chain lattice. In much the same 
way that tuning j'3 in the octahedral lattice produced helimag- 
netic states around the phase boundaries of the more common 
states, so in the chain lattice do we find that the asymmet- 
ric conic state becomes stabilized around what would be the 
antiferromagnetic, ferrimagnetic splayed, helimagnetic triple 
point. And unlike the helimagnetic state in the octahedral lat- 
tice, which had both first order and second order transitions, 
the asymmetric conic has only second order transitions (this is 
because it is an encompassing parametrization of every other 
state in this slice of the phase diagram). 

If we switch the sign of j 3 so that we have both ji = —1 
and j3 = —0.28 antiferromagnetic, the topology of the phase 
diagram (not shown) is much the same as Figure [8] The 
AFM/Heli boundary gets shorter and moves towards the lower 
left, so that the four domains almost meet at a point. More 
significantly, the asymmetric conic does not appear along the 
phase boundaries. Instead, all transitions are continuous, ex- 
cept that along both parts of FerriSplayed/Heli boundary, the 
portion closest to the center is first-order. [51] (The point 
where the nature of the transition switches from first-order to 
continuous is thus of tricritical type.) 

Finally, we consider the case of j\ ferromagnetic and j± 



FIG. 8: (COLOR ONLINE) Chain Lattice Phase Diagram for ji = 
— 1, j3=0.28, J4=0, j4=0. First state with asymmetric conic, and 
therefore first instance of non-coplanar ground states in the chain 
lattice. White circles denote helimagnetic phase, black circles 
anti-ferromagnetic, crosses splayed ferrimagnetic, and filled squares 
(blue online) asymmetric conic, the phases determined by numeri- 
cally optimizing the variational form of the energy for these cou- 
plings. 
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FIG. 9: (COLOR ONLINE) Chain Lattice Phase Diagram for ji = 
— 1, j3=0, J4 = —0.28, ji=0. Filled squares denote conic spiral 
phases, of which both kinds exist in this parameter slice: symmetric 
conic in the upper half (blue online) and asymmetric conic in the 
lower half (red online). This is our first instance of an alternating 
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antiferromagnetic (Figure[9]i. This slice of parameter space is 
quite interesting as both types of conies are present. Further- 
more, the alternating conic now fills a relatively large region 
of parameter space. This is likely due to its highly non-linear 
dependence on j'4 (as a function of its variational parameters). 

VIII. CONCLUSION AND DISCUSSION 

To conclude the paper, we first review our principal results, 
and then assess how much of what we learned is transferable 
to other lattices. 



Namely, when a ground state could not be written exactly, 
but emerges as an instability of a closely related state - e.g. 
a helimagnetic phase out of the ferromagnetic one - we can 
expand around the latter state and solve for the couplings at 
which it goes unstable. 

The insufficiency of finding regions with the correct Qlx 
for the purpose of stabilizing a targeted magnetic structures 
is an important caveat for the LT approach. There are some 
combinations of modes which cannot be stabilized by any set 
of pairwise interactions. That is, we can make them be op- 
timal modes, but only if they are part of a degenerate family 
that includes other optimal LT modes that are unrelated by 
symmetry. 



A. Summary 

The highlights of this paper include both concepts and 
methods, as well as results specific to the octahedral lattice, 
which seems relatively amenable to non-coplanar states. We 



pay special attention (Sec. VIII A3 1 to commonalities in the 
positioning of non-coplanar states in the phase diagram vis-a- 
vis neighboring phases. 

Our overall focus had a flavor of reverse engineering, in 
that we try to ask which couplings gave a certain phase (or 
which gave any non-coplanar phase) - of course in order to 
do that, one must also understand the forwards question (given 
the couplings, what is the phase). In that sense, our work is an 
example of a "materials by design" philosophy, whereby ma- 
terials are tailored - e.g. by adjusting their chemical content 
- to have a combination of interactions leading to a desired 
state. 



1. Methods 

Our basic recipe to determine the ground state of a non- 
Bravais lattice was a two-step process (Sec. First, an ap- 
proximate ground state configuration is generated through it- 
erative optimization (Sec. |HB| i of a lattice, starting from a 
random initial spin configuration. From this result, an ideal- 
ized spin configuration is created. The idealized formulation, 
if it has parameters undetermined by symmetry, is then used 
variationally optimize the Heisenberg Hamiltonian (Sec. |IIC| , 
yielding also the energy per site as a function of parameters. 
When this has been carried out for each candidate phase, a 
phase diagram (Sec. |VIl| can be generated. 

We refined the basic recipe further using three additional 
concepts or tricks. First, although the eigenmodes of the 
coupling matrix do not (in the non-Bravais case) automati- 
cally lead us to the ground state, they sometimes do work 
and are always a useful guide. Second, if the magnetic struc- 
ture is layered, the three-dimensional octahedral lattice can 
be mapped to various (non Bravais) one-dimensional lattices; 
such "chain" lattices are much more tractable than the cuboc- 
tahedral one, but they still support non-coplanar states. Third, 
we applied perturbation analysis to find second-order phase 
boundaries (especially in combination with variational opti- 
mization, but also in the Luttinger-Tisza Fourier analysis). 



2. Results 

We identified two categories of non-coplanar arrange- 
ments (which could be defined for more general lattices 
than the octahedral one). First, there are the commensu- 
rate three-Q states, exemplified here by the highly symmet- 



ric "cuboctabedral" spin state (Sec. IV 1, A second general 
class includes several varieties of incommensurate "conic spi- 
ral" (Sec. [V]). The commensurate three-Q states were found 
for both short and long range couplings (i.e. Jj's limited to 
no more than second-nearest neighbors or extending beyond 
second-nearest neighbors), while incommensurate conic spi- 
rals were possible for our lattice only when there are some 
Ji's beyond the second neighbors. 

The conic spirals came in "alternating" or "asymmetric" 
subclasses (stable in different parts of { Ji} parameter space) 
according to whether the spin components along the rotation 
axis were the same in each layer or alternated. This sug- 
gested that an additional subclass should exist, the "trans- 
versely modulated conic spiral", which is non-uniform in each 
layer. In addition, we came across a "double-twist" ground 
state of high (but not cubic) symmetry, which has some com- 
monalities with the cuboctahedral state but is probably best 
classified as a transverse-modulated conic spiral (Sec.|VI|>. 

The conic spiral states were found to be a function of one 
spatial coordinate, allowing the more involved octahedral lat- 
tice to be mapped to the ID chain lattice while preserving the 



ground state (Sec. V B 1. While this mapping is not a generic 
property of ground states in the octahedral lattice, it dramat- 
ically simplifies the analysis. Moreover, the transformation 
provides a guide for the couplings required to stabilize the 
conic spiral states in the octahedral lattice. 

Note that while the presence of inequivalent sublattices was 
necessary for the stabilization of non-coplanar states in the 
one-dimensional chain lattice, it is not so in higher dimen- 
sions. Our octahedral-lattice couplings explicitly treat all sites 
equivalently; non-coplanar states emerge either when it di- 
vides into unequal layers, that map to the chain lattice, or 
when it supports three-dimensional spin patterns not map- 
pable to a chain lattice. 

The most dramatic examples when the chain lattice fails 
to represent the ground state in the octahedral lattice are the 
cuboctahedral states, for which there is no distinguished di- 
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rection of variation (or stacking vector) 



1. Generalizing the cuboctahedral state? 



3. Why noncoplanar states are scarce 

One might have expected noncoplanar ground states to be 
generic (in the non-Bravais case, when they are possible at 
all) but in fact they were seen in only a small portion of our 



parameter space (Sec. VII I. In fact, non-coplanar phases typi- 
cally appeared in the phase diagram as "bridges" intermediate 
between simpler phases; there are two scenarios of bridging, 
as formulated in Sec. Ill Dl 

First, typically along the phase boundary between two sim- 
ple commensurate phases, one finds families of continuously 
degenerate states. Usually, this infinite family includes copla- 
nar states but the generic member is non-coplanar; appropriate 
tuning of the couplings (e.g. including further neighbors) can 
select particular non-coplanar states in different ways. How- 
ever, since the non-coplanar phase was limited (in the first- 
order description) to the (measure zero) boundaries between 
planar states, it naturally occupies only in a small region of 
the extended parameter space. Indeed, the non-coplanar states 
we found in the octahedral lattice either occurred in small 
wedges, e.g. the cuboctahedral or double-twist states or were 
only observed as accidental instances of the degenerate fam- 
ily in cases where the degeneracy cannot be broken (within 
the parameter space we took as the scope of this paper), e.g. 
the asymmetric conic states. 

Second, there are the "encompassing states". Such states 
become qualitatively different as a free parameter reaches 
some limiting value (a canonical example is the collinear fer- 
romagnet as a limiting case of planar spirals.) In the chain lat- 
tice, we often observed phase boundaries between two states 
that did not encompass each other; by tuning the couplings, it 
was sometimes possible to stabilize a third state which encom- 
passed the original two, so the original first-order transition is 
converted to two successive second-order transitions with the 
encompassing state in the middle. This is the means by which 
conic spiral states are stabilized in the chain lattice. 

Recently, a framework was proposed with a motivation sim- 
ilar to ours, the "regular states" ||52l which have a magnetic 
symmetry such that all sites symmetry-equivalent. The cuboc- 
tahedral state is an example. However, this is neither nec- 
essary nor sufficient for our own problem, to find all kinds 
of non-coplanar states: most of those, e.g. the double-twist 
state, are not regular, and conversely certain non-coplanar reg- 
ular states are not stabilized by bilinear exchange interactions 
alone. 



B. Generalizations to other lattices? 

Now we consider moving beyond the somewhat artificial 
octahedral lattice. We examine the necessary conditions for 
non-coplanar spin configurations analogous to the two main 
classes we discovered (cuboctahedral and conic spiral). Can 
our results be applied to other lattices, such as the pyrochlore 
lattice? 



The cuboctahedral states seem highly specific to the octa- 
hedral lattice, since they possess the same symmetries as the 
lattice. On the other hand, the LT construction (each sublattice 
using different combinations of LT modes, the combinations 
being related by rotational symmetry of the different sublat- 
tices) seems fairly generic. A general name for states like the 
cuboctahedral state might be "commensurate triple-Q state", 
referring to its content of LT modes. Can such states be found 
in other lattices, or are there others in the octahedral lattice? 

One answer requires that the hypothetical generalized state 
enjoys the full symmetry of the lattice, as the cuboctahedral 
state does. Then the number of site classes must be a multiple 
of the number of spin components, i.e. of three. That will not 
work for the pyrochlore lattice, in which there are four site 
classes associated with (111) directions. However, it might 
work in the half-garnet lattice, which has six site classes asso- 
ciated with (1 10) directions. 

A second, more systematic way to answer question follows 



the LT approach of Sec. IV C The state must be a linear com- 
bination of three modes - one for each component of spin. For 
the state to enjoy the full lattice symmetry, these must be the 
complete star of symmetry-related modes, thus it must have a 
threefold multiplicity. Furthermore, each mode must be real 
(otherwise one mode requires two spin components); that hap- 
pens only if the wavevector is half a reciprocal lattice vector, 
i.e. is at the center of one of the Brillouin zone's faces. 



However, Sec. IV C 1 gave a cautionary example: finding 
such modes is not sufficient, because it might give a state with 
decoupled sublattices, so that the cuboctahedral state would 
merely be one undistinguished configuration in a continu- 
ous manifold of degenerate states that even includes collinear 
states. We can, in fact, borrow a notion from Section fVBl to 
guess when this happens: mapping to a chain lattice, now ap- 
plied to LT modes rather than to spin structures. The wavevec- 
tor of one mode defines a set of planes, and thus a way of 
projecting both the sites and the mode onto a one-dimensional 
chain lattice. On the chain lattice, in order that the mode be 
real, the wavevector must be it. If all the site planes were 
equivalent, the chain lattice is a Bravais lattice and the mode 
must be a plane wave; with Q = tt, it is easy to see we get 
a decoupled pattern, namely (+1, 0, —1, 0, ...). In the octa- 
hedral lattice, the {1/2, 0, 0} modes have a threefold degen- 
eracy and also correspond to a non-equivalent set of stacked 
planes, giving a cuboctahedral state, whereas the threefold 
degenerate {1/2,0,0} modes have equivalent stacked planes 
and give decoupled states. In the pyrochlore lattice, the 
{1, 0, 0} modes have a threefold multiplicity but have equiv- 
alent stacked planes; on the other hand, the {1/2, 1/2, 1/2} 
modes have inequivalent stacked planes but their multiplicity 
is fourfold. Thus, in the pyrochlore case this path does not 
lead us to a cuboctahedral state. 

The third way to answer this question is via the cluster con- 



stuction of Sec. IV A in which the lattice was decomposed 
into cuboctahedral cages. Indeed, some other lattices consist 
of a union of roughly spherical "cage" clusters: e.g., in the 
CraSi-type or "A15" structure, the majority atoms form cages 
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in the form of distorted icosahedra. Furthermore, couplings 
can certainly be chosen such that the spin ground state of a 
cluster forms a the same polyhedron as the cluster itself. 



However, our ground-state construction (Sec. IV A i de- 
manded that the spin configurations in adjacent cages be re- 
lated by a reflection in spin space. Thus the cages alternate 
between proper and improper rotations of a reference config- 
uration. That is possible only if the cage centers form a bi- 
partite lattice, which is not true for the A15 majority-atom 
cage centers. (They form a bcc lattice with first- and second- 
neighbor links.) 



2. Generalizing the non-coplanar spirals ? 

Long ago, Kaplan and collaborators studied doubly mag- 
netic spinels, and identified a (noncoplanar) "ferrimagnetic 
spiral" configuration arising from exchange interactions (a 
spiral of this sort is responsible for multiferroic properties of 
CoCr2C>4 lfl8l . This is actually a kind of double spiral (each 
magnetic species accounting for one of the spirals). Notice 
that we could map that structure to a chain, too, and that chain 
would have two inequivalent sites, which we point out is a 
precondition for developing a noncoplanar state. 

For a mapping to the ID chain lattice to provide a non- 
coplanar state, the mapped sites must be inequivalent by trans- 
lation; if not, the chain lattice is a Bravais lattice and must 
have (at most) a coplanar spiral state. This in turn depends on 
having unequal layers in the three-dimensional lattice, which 
is is impossible in the Bravais case, but inequivalent layers can 
emerge from a non-Bravais lattice wit h full y symmetric cou- 
plings. Thus, as we worked out in Sec. V B a {100} stacking 
in the octahedral lattice has twice as many sites on x and y 
bonds, constituting one kind of layer, as there are sites on z 
bonds, constituting the other kind. 

The stacking direction is one in which the wavevector is 
incommensurate, and so we may seek out parameter sets 
for which Qlt goes incommensurate in a desired direction. 
However, this is no guarantee that the actual ground state is 
stacked, since it might be better energetically to combine these 
modes (or ones nearby in k-space) in a quite different way. 
Thus, although we found some optimal Qlt wavevectors in 
directions other than {100}, they were never the basis of a 
stacked spiral. 



C. On to pyrochlore lattice? 

The pyrochlore lattice is undoubtedly the most prominent 
non-Bravais lattice that is likely to have non-coplanar states, 
and was a major motivation for the methods we developed in 
this paper for the more tractable octahedral lattice. 

Relatively few works have tackled the isotropic spin prob- 
lem on the pyrochlore lattice, beyond the (massively degen- 
erate) case of only J\ [analogous to the J\ = J2 case on the 
cuboctahedral lattice]. The J1-J2, or J1-J3 pyrochlore, with 
large antiferromagnetic Jj, has a noncoplanar and somewhat 
obscure state that is not fully understood at T = 



In the spinel GeNi204, neutron diffraction found a {111} or- 
dering consisting of alternating kagome and triangular lattice 
layers with different densities of spins, reminiscent analogous 
of the {100} stackings we have explored in the octahedral lat- 
tice. Ref. [53] invoked interactions out to J4 in order to ra- 
tionalize this state, but did not verify it was the ground state 
for the suggested interactions. Finally, to address metallic py- 
rochlore compounds, a toy Hamiltonian was investigated on 
the pyrochlore lattice with exchange couplings having RKKY 
oscillations l57l . 

Very recently, a comprehensive study has been carried out 
for just the J\-Ji phase diagram of the pyrochlore lattice by 
applying the methods of the present paper, which were equally 
effective in that case ||58ll . About four complex states were 
identified, including that above mentioned state introduced by 
Kawamura and collaborators [ 54-56 1; a state reminiscent of 
our double-twist state; a more complex generalization of a 
conic spiral; and a cuboctahedral spin state that has less than 
cubic lattice symmetry, thus evading the negative conclusions 
of Sec.[VTinn]|58). 



Acknowledgments 

We thank G.-W. Chern, J-C. Domenge, C. Fennie, S. 
Isakov, T. Kaplan, M. Lapa, and M. Mostovoy for discus- 
sions and communications. C.L.H. was supported by NSF 
grant DMR- 1005466; S. S. was supported by the Cornell En- 
gineering Learning Initiatives Program, the Intel Corporation, 
and the Rawlings Cornell Presidential Research Scholars Pro- 
gram. 



Appendix A: Perturbative Calculation of Phase Boundaries 

To move from a collection of ground states discovered at 
discrete points in parameter space to draw a full phase dia- 
gram is non-trivial. To analyze the phase boundaries of the 
more complicated states (those with free parameters in their 
parametrization) we depend on either variational optimization 
(Sec. refsec:var), or some kind of perturbation theory. Per- 
turbation theory can be applied in two places: either to the 
LT matrix, or to the Hamiltonian of a parametrized spin state. 
The former is more straightforward, but is limited since most 
of the non-coplanar ground states are not exactly built from 
optimal LT eigenmodes. 

An obvious caveat for either of these applications of pertur- 
bation is that they detect continuous transitions, representing 
infinitesimal changes in the spins: it is a bifurcation of the lo- 
cal minima as points in the configuration space. But there the 
ground state might instead change due to a first order transi- 
tion, when the energies of two separated configurations cross 
as parameters are varied. We do encounter the first-order case 
on occasion, though not nearly so often as the continuous one. 
To detect such discontinuous transitions, we must compare 
numerical calculations of the ground state energies. 

Our original question was "given a certain set of param- 
eters, what is the ground state", but in these calculations it 
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has been turned around to "given a particular ground state, for 
what parameter sets is it favored?" 



1. LT Analysis (octahedral lattice example) 

We give an example here of the use of perturbation theory 
to discover the incipient instabilities of LT modes. Such an 
approach may be especially useful to locate the phase bound- 
aries for transitions from an LT exact ground state to more 
complicated, incommensurate states. If one is hunting for the 
parameter domain which would stabilize an particular mode 
with ordering wavevector Q, a substitute problem is to find 
the parameter domain in which this mode (or one with Qlt 
similar to Q) is the optimal LT mode. We reiterate the caveat 



from Sec. II A such a discovery is a necessary but not suffi- 
cient criterion to guarantee the existence of any ground state 
based on the obtained Qlt- (It is not sufficient because the 
actual ground state could feature additional modes or modes 
merely in the neighborhood of Qlt-) 

The LT matrix elements Jy(k) and eigenvalues J(kv) are 
functions of wavevector k. Imagine that Qo is a point of high 
symmetry in the zone so as to be a stationary point for J(kz^). 
For some parameter sets, we know, it is a minimum and in fact 
optimal; whereas for some other parameter sets, we imagine, 
it is only a saddle point, and the minimum occurs at some 
nearby wavevector of lower symmetry. 

One first writes a Taylor expansion of the LT matrix in pow- 
ers of (5k = k Qo- Using standard techniques (formally 
identical to those used for eigenfunctions of the Schrodinger 



equation) it is straightforward to write a perturbation expan- 
sion for J(ki/) in powers of <5k. Inspection then shows where 
this stops being positive definite. Since the LT matrix ele- 
ments are bounded, so are the eigenvalues J(ki/). So if the 
mode at Qo goes unstable at quadratic order in (5k, there must 
be higher-order positive terms in the expansion. Thus the sin- 
gle local minimum of J{kv) bifurcates in some fashion. The 
corresponding spin state cannot be a commensurate spiral, but 
it might be representable in the framework of planar stackings 
(Sec.|VB|>. 

This technique was used, for example, to analyze how the 
Qo = (1/2, 1/2, 0) wavevector is destabilized in the LT phase 
diagram. Simulations had found degenerate antiferromagnetic 
orderings at that wavevector; if an incommensurate wavevec- 
tor of form (q, q, 0) had been stabilized, this might have been 
the basis of a non-coplanar spiral stacked in the (1,1,0) di- 
rection. 

The LT matrix is given by: 



J«( k ) 



J y (k) 



2 J2 cos ki + 2 J^cos kj + cos kk) + 4 J 4 cos ki 

x (cos kj + cos kk) + cos kj cos kk 

(Ala) 

k ' k ~ k ' k ~ 

4Ji cos — cos + 8 J3 cos -j- cos cos kk- 

(Alb) 



Here j, k in (Alai means the indices other than i, similarly 
k in ( |Alb| i is the index other than i, j. If we substitute k = 
(k, 0, 0), for example, the eigenvalues along this cut are 



J(k, 0) = - J 2 + 24 - (J F (0) - J 2 )[l + cosfc] + 2J x (k); 
J(k,±) = --[24-4J 4 + (Jir(0)+4J 4 )(l + cosfe) + 2Jx(fc) 



± 2 \ K J 2 - J 2 + 2J 4 - 24)(1 - cosfc) + 2J^(fc)] 2 + 16J|(0)(1 + cosfc) 



1/2 



(A2a) 



(A2b) 



Along another cut through the zone, k = (fc, fc, 0)/2y / 2, the eigenvalues are 



J(k, 0) = J 2 - 2J' 4 - 2J F {0) cos 2 fc + 2J X (0) cos 2 fc; 

J(k,±) = J' 2 + 2J 4 cos 2 2fc - (Jp(2fc) + 2J' 2 + 4 J 4 cos 2fc) cos 2 fc - 2J x (0)cos 2 fc 
± 2^ [(*M2fc) - 4 J 4 cos 2fc) sin 2 fc - J x (0)cos 2 fc] 2 + 8J x (2k) 2 cos 2 fc) 



1/2 



(A3a) 



(A3b) 



In both cuts, where Jf(?) = Ji + J'2 + 2 J 4 cos q + 2 J4 cos q is the effective ferromagnetic coupling and Jx (q) = Ji + 
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2J3 cosq is the effective cross-sublattice coupling. 



2. State Perturbation (chain lattice example) 

The major advantage of applying perturbation theory to a 
parametrization of the lattice's spins, rather than the LT ma- 
trix, is that it can accomodate a more generic ground state. 
This is not say that it avoids the disconnect between the state 
being perturbed and the actual ground state, if the perturbed 
state does not encompass the true ground state then the dis- 
connect is required. However, this method does allow us to 
consider states composed of multiple LT wave-vectors. 

In the octahedral lattice, most of the ground state types we 
encountered - even the noncollinear ones - are essentially 
built using a single "star" of symmetry-related ordering wave- 
vectors. Therefore, for the octahedral lattice, the actual phase 
diagram mostly reflects the LT phase diagram and it is prefer- 
able to find phase boundaries via the LT perturbation method 
illustrated above in Sec. |A 1| For the chain lattice, however, 
the ground state is typically characterized by several wave- 
vectors. 

(as noted near the beginning of Sec. [V]) 

Such multi-LT wave-vectors are all parametrized by some 
form of "conic spirals," a mix of a planar spiral using the 
wavevector (q) and a deviation along the wavevector (£2) 
where fc 2 is either integer or half-integer (since the chain lat- 
tice is ID, the wave-vectors are as well). Within a range of 
parameter space j\ through (75 an d higher all 0), these 
conic spirals are the most general form of ground state, mak- 
ing the problem of finding the ground state amenable to vari- 
ational methods. 

As an illustration of determining the ground state by varia- 
tional methods, we consider the problem of finding the phase 
boundaries for the "alternating conic" k 2 — 1/2. The spin 



configurations in the lattice are parametrized by Eq.(5.3bi. 



Without loss of generality, we take the odd sites to be the ones 



with planar spins. A symmetry relates the states (5.3b 1 to the 
other family of alternating conic configurations in which even 
and odd sites have swapped roles, if one also swaps in param- 
eter space (j2,ji) «-> (JqJ'a)- 

Up to interaction j 4 , the energy per unit cell is 

E = —2{ji cos i> + j 3 cos 3V>) cos a 

— 2(j2 cos 2 ib — H sin 2 20) cos 2 a 

+ ]2~ ji- i' 2 cos2il> - j' 4 cosAip (A4) 



Setting dE/d(cosa) = to minimize (A4i, we see that the 
optimal angle a*, is given by 



cos a = — - 



1 ji+j 3 (4cos 2 y>-3) 
2 cos V> J2 - 4j 4 sin 2 if) 



(A5) 



Notice the symmetry under reversing the signs of ji and j'3 
and a f-> tt — a. Of course, a necessary condition is that the 
r.h.s. of ((AS} lies in (0, +1) (recall that < a < ir/2 by 



definition), otherwise a* is pinned to or tt/2, which would 



be a planar spiral (or possibly a collinear state, depending on 
the value of ip). 



Plugging (A5i into Eq. (A4i leaves ip = q/2 as the only 
variational parameter: 



E = 



J ^=^ + cos 2 ib 



j 2 cos 2i/j 



-j' 4 cos 40 - 32 - h 



(A6) 



If we drop the secondary couplings equation (A 2 1 reduces to 

j'2 cos 2 V> 



E 



2 32 



(A7) 



Since ( A7 1 is linear in cos 20, its ground state given by 
ip m = mn/2, with m integer. If m is odd, then Eq |A7| breaks 
down, since Eq.( A5 1 is singular. What actually occurs here is 
that the sublattices have adopted an antiferromagnetic struc- 
ture, decoupling the sublattices. Because the sublattices are 
decoupled, a is arbitrary (representing the freedom the de- 
coupled sublattices to rotate relative to each other). The en- 
ergy for such a configuration is E = 32 + j' 2 - Conversely, if 
m is even, then the ground state is a splayed state (assuming 
that < \2j 2 \) That means a commensurate, planar state 
(using the xz plane, or collinear if a* is trivial). Thus, with 
only primary couplings, even though the two sublattices are 
inequivalent, we cannot obtain a non-coplanar spiral. 

When we turn on (not too large) secondary couplings, those 
commensurate states will remain stable out to the critical cou- 
pling, at which the optimal wavevector bifurcates. Therefore 
we expand ( A 2 1 in powers of 8 = tp — ip m about both kinds 



of stationary point [59] to find: 

j^,^ (h + J3) 2 ./ ./ 

E{o) = — j 2 - Ji - 32 - J4+ 



2j2 

ro 0i +J3) 2 . , r . 



J2 



-8j 4 ]<5 2 
(A8a) 



V 2 ' 2 j 2 -4i4 



+ fa - 34 ~ 32 - 34+ 



. 3i - 3j 



32 - 4j 



-.73 + 2[ 



3i - 3.73,2 



32 - 4j 



-} 2 34 - 23 / 2 + S3 / i \S 2 



(A8b) 



The commensurate state becomes unstable to S ^ when the 
coefficient of the quadratic term goes negative, so the condi- 
tions to induce instability are: 

E{d):0> ^ 34-2 J3+J 2 + 4 J4 



32 



(A9a) 



Ek : > 



,3i - 3j 



3l2, 



32 ~ 4 J4 



34 



, Jl - 3j; 

' 32 - 4jz 



.73 - 32 + 4 J4 



(A9b) 
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Combined with the requirement that |cosa| < 1 in (A5i, istance of an alternating conic state. 
Eqs. (|A9]> give the minimum necessary conditions for the ex- 



[1] A.Yoshimori, J. Phys. Soc. Jpn. 14, 807 (1959); J. Villain, J. 

Phys. Chem. Solids 11, 303 (1959); T. A. Kaplan, Phys. Rev. 

116, 888 (1959). 
[2] T. A. Kaplan and N. Menyuk, Phil. Mag. 87, 371 1 (2006) 
[3] C. L. Henley, Phys. Rev. Lett. 62, 2056 (1989). 
[4] E. F. Shender, Sov. Phys. JETP, 56, 178 (1982). 
[5] C. L. Henley and B. E. Larson, preprint (arxiv 2008); B. E. 

Larson and C. L. Henley, "Ground State Selection in Type III 

FCC Vector Antiferromagnets" (unpublished preprint, 1990). 
[6] X. G. Wen, F. Wilczek, and A. Zee, Phys. Rev. B 39, 11413 

(1989). 

[7] S. Sachdev (personal communication). 

[8] L. Messio, B. Bernu, and C. Lhuillier, Phys. Rev. Lett. 108, 
207204 (2012). 

[9] Y. Taguchi, Y. Oohara, H. Yoshizawa, N. Nagaosa, and Y. 
Tokura, Science 291, 2573 (2001), 
[10] A. Kalitsov, B. Canals, and C. Lacroix, J. Phys. Conf. proc. 145, 
012020 (2009). 

[11] M. Taillefumier, B. Canals, C. Lacroix, V. K. Dugaev, and P. 

Bruno, Phys. Rev. B 74, 085 105 (2006) 
[12] N. Nagaosa, J. Sinova, S. Onoda, A. H. MacDonald, and 

N. P. Ong, Rev. Mod. Phys. 82, 1539 (2010) 
[13] C. L. Henley, Ann. Phys. (N. Y.) 156, 368-411 (1984), 
[14] S.-W. Cheong and M. Mostovoy, Nature Materials 8, 13 (2007). 
[15] T. Kimura, Ann. Rev. Mater. Sci. 37, 387 (2007). 
[16] T. Kimura, J. C. Lashley, and A. P. Ramirez, Phys. Rev. B 73, 

220401 (2006) 
[17] D. Khomskii, Physics 2, 20 (2009). 

[18] Y. J. Choi, J. Okamoto, D. J. Huang, et al. K. S. Chao, H. J. Lin, 
C. T. Chen, M. van Veenendaal, T. A. Kaplan, and S. W. Cheong 
Phys. Rev. Lett. 102, 067601 (2009). 

[19] R. Moessner and J. T. Chalker, Phys. Rev. Lett. 80, 2929 (1998); 
Phys. Rev. B 58, 12049 (1998). 

[20] S. T. Chui, Phys. Rev. B 15, 307 (1977). 

[21] P. Reed, J. Phys. A 10, 1745 (1977). 

[22] M. Hermele, M. P. A. Fisher, and L. Balents, Phys. Rev. B 69, 
064404 (2004). 

[23] T. S. Pickles, T. E. Saunders, and J. T. Chalker, EPL 84, 36002 
(2008). 

[24] R. Shankar, F. J. Burnell, and S. L. Sondhi, Ann. Phys. (N.Y.) 
324, 267 (2008). 

[25] H. Takizawa, T. Yamashita, K. Uheda, and T. Endo, J. Phys. 

Condens. Matter 14, 1 1 147 (2002) 
[26] D. Fruchart, E. F. Bertaut, J. P. Senateur, and R. Fruchart, J. de 

Physique Lett. 38, L-21 (1977); D. Fruchart and E. F. Bertaut, 

J. Phys. Soc. Jpn. 44, 781 (1978) 
[27] T. Koyama, H. Yamashita, Y. Takahashi, T. Kohara, I. Watan- 

abe, Y. Tabata, and H. Nakamura, Phys. Rev. Lett 101, 126404 

(2008). 

[28] V. H. Tran, W. Miller, and Z. Bukowski, Phys. Rev. Lett. 100, 
137004 (2008). 

[29] T. Nagamiya, "Helical Spin Ordering: 1 Theory of Helical Spin 
Configurations", in Solid State Physics vol. 20, p. 306, ed. F. 
Seitz, D. Turnbull, and H. Ehrenreich (Academic Press, New 
York and London 1967). 

[30] J. M. Luttinger and L. Tisza, Phys. Rev. 70, 954 (1946). 

[31] D. H. Lyons and T. A. Kaplan, Phys. Rev. 120, 1580 (1960). 



[32] J. Samuel Smart, Effective field theories of magnetism (W. B. 
Saunders, Philadelphia, 1966); see chapter 8. 

[33] D. Bergman, J. Alicea, E. Gull, S. Trebst, and L. Balents, Na- 
ture Physics 3, 487 (2007). 

[34] S. Okumura, H. Kawamura, T. Okubo, Y. Motome, J. Phys. Soc. 
Jpn. 79, 1 14705 (2010); A. Mulder, R. Ganesh, L. Capriotti, and 
A. Paramekanti, Phys. Rev. B 81, 214419 (2010). 

[35] D. H. Lyons, K. Dwight, T. A. Kaplan, and N. Menyuk, Phys. 
Rev. 126, 540 (1962). 

[36] M. J. Freiser, Phys. Rev. 123, 2003 (1961). 

[37] L. R. Walker and R. E. Walstedt, Phys. Rev. Lett. 38, 514 (1977) 

[38] C. L. Henley, Ann. Phys. (N. Y.) 156, 324 (1984). 

[39] C. L. Henley, Can. J. Phys. 79, 1307 (2001). 

[40] V. Elser, I. Rankenburg, and P. Thibault, PNAS 104, 418 (2007). 

[41] This is not quite the magnetic structure factor, since that in- 
cludes interferences between different sublattices. 

[42] Occasionally the end-to-end plot is illuminating even in higher 
dimensions, the sequence of spins being defined by a raster scan 
that goes down rows - limited to one of the three sublattices 
- running in one of the coordinate directions.) For example, 
the double-twist state (Sec, |VIfr was discovered from common- 
origin and end-to-end plots. 

[43] E. Rastelli, A. Tassi, and L. Reatto, Physical B & C 97, 1 
(1979); E. Rastelli and A. Tassi, J. Phys. C: 19, L423 (1986). 

[44] J. N. Reimers, A. J. Berlinsky, and A. C. Shi, Phys. Rev. B 43, 
865 (1991). 

[45] D. H. Lyons and T. A. Kaplan, J. Phys. Chem. Solids, 25, 645 
(1964). 

H.-J. Schmidt and M. Luban, J. Phys. A 36, 6351 (2003); 

[46] J.-C. Domenge, P. Sindzingre, C. Lhuillier, and L. Pierre, Phys. 
Rev. B 72, 024433 (2005). 

[47] J.-C. Domenge, C. Lhuillier, L. Messio, L. Pierre, and P. Viot, 
Phys. Rev. B 77, 172413 (2008). 

[48] Spin configurations on a cuboctahedron were studied by: J. 
Schulenberg, A. Honecker, J. Schnack, J. Richter, and H.-J. 
Schmidt, Phys. Rev. Lett. 88, 167207 (2002). 

[49] We did not include the asymmetric conic state in table [I] since 
we will see that it is only found on the phase boundary between 
two other states. For the octahedral lattice, then, it appears only 
as a randomly chosen member of a family of highly degener- 
ate states (using the classification of Sec. |IID} , rather than a 
genuine ground state. 

[50] Due to the complexity of the double-twist state, the phase 
boundaries of the double-twist state shown in Figure|6]were not 
generated by variationally optimizing a parametrized Hamilto- 
nian and finding absolute minima. Instead, we show the phase 
boundary from the "Luttinger-Tisza phase diagram", i.e. the 
boundary of the region of parameter space in which (q, q, 0) 
was the optimal wave-vector. As argued in Sec. |II A| that will 
produce a phase diagram topologically similar to the true phase 
diagram. 

[51] Around the tricritical-type points where the transition switches 
from first-order to degenerate, there are alternating conic states 
only slightly higher in energy than the helimagnetic states, so 
that a sliver of alternating conic may be spuriously identified 
here if numerical results are not checked analytically. 

[52] L. Messio, C. Lhuillier, and G. Misguich, Phys. Rev. B 83, 



23 



184401 (2011). 

[53] M. Matsuda, J.-H. Chung, S. Park, T. J. Sato, K. Matsuno, H. 

Aruga-Katori, H. Takagi, K. Kakurai, K. Kamazawa, Y. Tsun- 

oda, I. Kagomiya, C. L. Henley, and S.-H. Lee, Europhys. Lett. 

82, 37006 (2008). 
[54] D. Tsuneishi, M. Ioki, and H. Kawamura, J. Phys. Condens. 

Matter 19, 145273 (2007). 
[55] T. Nakamura and D. Hirashima, J. Mag. Mag. Materials 310, 

1297 (2007). 

[56] G.-W. Chern, R. Moessner, and O. Tchernyshyov Phys. Rev. B 



78, 144418 (2008) . 

[57] "Ordering of the pyrochlore Ising model with the long-range 
RKKY interaction" A. Ikeda and H. Kawamura J. Phys. Soc. 
Jpn. 77, 073707 (2008). 

[58] M. Lapa, B.S. thesis, Cornell University (2012); M. Lapa and 
C. L. Henley, in preparation. 

[59] There are just two kinds of stationary point since ( |A7| > has pe- 
riod 7T. 



